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Abstract 

We study spectral asymptotics for small non-selfadjoint perturbations of 
selfadjoint /i-pseudodifferential operators in dimension 2, assuming that the 
classical flow of the unperturbed part possesses several invariant Lagrangian 
tori enjoying a Diophantine property. We get complete asymptotic expansions 
for all eigenvalues in certain rectangles in the complex plane in two different 
cases: in the first case, we assume that the strength e of the perturbation is 
0(h s ) for some 5 > and is bounded from below by a fixed positive power of 
h. In the second case, e is assumed to be sufficiently small but independent 
of h, and we describe the eigenvalues completely in a fixed /i-independent 
domain in the complex spectral plane. 
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1 Introduction and statement of main results 

Recently there has been a large number of new developments for non-selfadjoint 
problems. These include semiclassical spectral asymptotics for non-selfadjoint ope- 
rators in low dimensions [22], 121], [21] • [10], [US], resolvent estimates and pseudospec- 
tral behavior [TH], [IB], M spectral instability questions [TO], [SI], and evolution 
problems and decay to equilibrium for the Fokker-Planck operator [20J. The pur- 
pose of this work is to continue a line of development initiated in |28j . which opened 
up the possibility of carrying out a spectral analysis for non-selfadjoint operators in 
dimension two, that is as precise as the corresponding analysis for selfadjoint prob- 
lems in dimension one. In [28 , it was established that for a wide and stable class of 
non-selfadjoint operators in dimension two, it is possible to describe all eigenvalues 
individually in a fixed domain in the complex plane, by means of a Bohr- Sommerf eld 
quantization condition. The underlying reason for this result is a version of the KAM 
theorem without small divisors, in a complex domain. 

The work [2E] has been continued in a series of papers [Jl] , [22] , [23] , [2H] , all of 
them done in the context of small non-selfadjoint perturbations of selfadjoint opera- 
tors, with the important additional assumption that the classical flow of the leading 
symbol of the unperturbed part should be periodic in some energy shell. While the 
case of a periodic classical flow is very special indeed, in the aforementioned works, 
we have already given some applications of the general results to spectral asymp- 
totics for damped wave equations on analytic Zoll surfaces [30], [21], while barrier 
top resonances for semiclassical Schrodinger operators have been treated in |2*o] . 

Now a classical Hamiltonian with a periodic flow can be naturally viewed as a 
degenerate case of a completely integrable symbol, and an even more general and 
much more interesting dynamical situation occurs when considering a symbol that 
is merely close to a completely integrable one. Continuing our previous works, in 
this case it seems to be of interest to study the spectrum of non-selfadjoint operators 
that are small perturbations of a selfadjoint operator, whose classical flow is close 
to a completely integrable one. The present work is the first one where we begin 
to study this problem, and when doing so, as our starting point, we shall take a 
general assumption that the real energy surface of the unperturbed leading symbol 
contains several flow invariant Lagrangian tori satisfying a Diophantine condition. 
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According to a classical theorem of Kolmogorov |3] , the existence of such invariant 
tori is guaranteed when the unperturbed symbol in question is close to a completely 
integrable non-degenerate one. 

We shall begin by describing the general assumptions on our operators, which 
will be the same as in (221 > EI|> an d |2H]- Let M denote R 2 or a compact real- 
analytic manifold of dimension 2. We shall let M stand for a complexification of M, 
so that M = C in the Euclidean case, and in the manifold case, M is a Grauert 
tube of M. 

When M = R 2 , let 

P e = P(x,hD x ,e;h) (1.1) 

be the Weyl quantization on R 2 of a symbol P(x, £, e;h) depending smoothly on 
e G neigh (0, R) with values in the space of holomorphic functions of (x,£) in a 
tubular neighborhood of R 4 in C 4 , with 

\P(x,Z,e;h)\<O(l)m(Re(x,0) (1.2) 

there. Here m is assumed to be an order function on R 4 , in the sense that m > 
and 

m(X) <C (X- Y) N "m(Y), X, Y e R 4 , C , N > 0. (1.3) 
We also assume that 

m>l. (1.4) 

We further assume that 

oo 

P(x,£,e;/i)~^- e (:r,0^, h ^ 0, (1.5) 
i=o 

in the space of such functions. We make the ellipticity assumption 

K,(x,OI > ^m(Re(x,0), \(x,£)\>C, (1.6) 

for some C > 0. 

When M is a compact manifold, we let P e be a differential operator on M, such 
that for every choice of local coordinates, centered at some point of M, it takes the 
form 

Pe=J2 a<*,e( x ;h)(hD x ) a , (1.7) 

|a|<m 
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where a a ^(x;h) is a smooth function of e with values in the space of bounded 
holomorphic functions in a complex neighborhood of x — 0. We further assume 
that 

oo 

a ate (x;h) ^ J~]g ate j(x)h 3 , h -»• 0, (1.8) 
j"=o 

in the space of such functions. The semi-classical principal symbol po,e, defined on 
T*M, takes the form 

P0,e( X i = Y1 a a,efi( X )C (1-9) 

if (x, £) are canonical coordinates on T*M, and we make the ellipticity assumption 

bo,e(^OI > ^<0 m , foO e lei > c, (1.10) 

for some large C > 0. (Here we assume that M has been equipped with some 
Riemannian metric, so that |£| and (£) = (1 + l^] 2 ) 1 / 2 are well-defined.) 
Sometimes, we write p e for p t and simply p for p . Assume 

P e=0 is formally selfadjoint. (1-H) 

In the case when M is compact, we let the underlying Hilbert space be L 2 (M, fx(dx)) 
for some positive real-analytic density /i(dx) on M. 

Under these assumptions, P e will have discrete spectrum in some fixed neighbor- 
hood of G C, when h > 0, e > are sufficiently small, and the spectrum in this 
region will be contained in a band \lmz\ < 0(e). 

Assume for simplicity that (with p = p t= o) 

p-\0) n T*M is connected, (1.12) 

and let us also assume that the energy level E = is non-critical, so that dp ^ 
along p _1 (0)nrM. 

Let Hp = p^ ■ — p' x ■ ^ be the Hamilton field of p. We introduce the following 
hypothesis, assumed to hold throughout this work: 

The setp _1 (0) fl T*Mcontains finitely many analytic H p - invariant (1.13) 
Lagrangian tori Aj, 1 < j < L, such that each Aj carries real analytic 
coordinates x±, x 2 identifying Aj with T 2 so that along Aj,we have, 

H p = a ± d xi + a 2 d X2 , (1.14) 
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where a\, a 2 G R satisfy the Diophantine condition, 



\a-k\>—}-^, O^keZ 2 , (1.15) 

Gn \k\ 



for some fixed C , N > 0. Here T 2 = R 2 /2ttZ 2 is the standard 2-torus. 

We write out the first few terms in a Taylor expansion of p e in a neighborhood 
of p _1 (0) H T*M, 

p e =p + i eg + C'(e 2 ). (1.16) 

When < A" G C£° (R) is such that / K(t) dt = 1 and T > 0, we introduce a 
"smoothed out" flow average of q, 



(q)r,K,p = (q) 



T,K 



j K T (-t)q o exp (tH p ) dt, K T {t) = (lj , (1.17) 



defined near p _1 (0) H T*M. The standard flow average corresponds to taking K = 
l[-i,o], an d we shall then write {q)r,K = {q)t- 

Let Gt be an analytic function, defined in a neighborhood of p _1 (0) fl T*M, such 
that 

H P G T = q~ (q)r,K- 

This is a convolution equation along the if p -trajectories, and as in jSHj and [21], we 
solve it by setting 



(1.18) 



Gt = j TJ T {-t)q o exp {tH p ) dt, J T (t) = i J , 

where the function J is compactly supported, smooth away from 0, with 

J'{t) =5(t)-K(t). 

Composing the principal symbol ()1.16j) with the holomorphic canonical transfor- 
mation exp (ieHa T ), and conjugating the operator P t by means of the corresponding 
Fourier integral operator U e = e^ GT ^ x ' hDx \ defined microlocally near p _1 (0) fl T*M, 
we may reduce our operator to a new one, still denoted by P e , which has the principal 
symbol 

p e o exp (ieH Gr ) =p + ie(q) TjK + T (e 2 ). 

Moreover, it is still true that P e =o is the original selfadjoint operator. Repeating 
an argument, explained for example in [HH!, which makes use of the sharp Garding 
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inequality, we obtain a first localization of the spectrum of P t : if z G C in the 
spectrum of P e is such that |Rez| < 5, then as e, S, h — > 0, 



Imz 
G 



lim inf Re (q)TK ~ nm SU P R- e (q)tk + o(l) 

T^oop-i(O) ' T ^°°p-i(0) 



;i.i9) 



This estimate remains valid for K = l[-i,o]- Let us also notice that along the 
diophantine torus Aj, 1 < j < L, we have uniformly, as T — > oo, 

<<z>t = + (~) • (1-20) 

Here Fj is the mean value of g over Aj , computed with respect to the natural smooth 
flow-invariant measure on Aj, with respect to which the if p -flow on Aj is ergodic. 
In the case when K G Cq°(R), using the rapid decay of K, it is easy to see that 
(jl.20|) improves to 

1 



(q) T , K = F J + ^— 

We shall assume from now on that 

F\ = F 2 = . . . = F L is independent of j, (1-21) 

and we write then Fj = F,l<j<L. 

As will be explained in section 2, for each j, there exists a smooth canonical 
transformation 

KooJ : neigh (Aj, T* M) -> neigh(£ = 0,T*T 2 ), 

mapping Aj to £ = 0, such that 

Pe o k"^ = PooM + m(x,Z) + 0(e 2 ) + 0(r), (1-22) 

and Poo,j(0 = a ' £ + C*(£ 2 )) with a satisfying ()1.15|) . We furthermore may assume 
that the energy surface p^j(0) has the form £ 2 = fj(£,i), for some smooth function 
fj, with /j(0) = 0, /j(0) 7^ 0. Using the coordinates £i, £ 2 > we define 



(^•)(0 = ^p y <^,o^- (i-23) 

For each small 5 > 0, we use the coordinate functions £i o k^j and £ 2 K oo,j near 
Aj to decompose the real energy surface as follows, 

p-^O) = fi_(<5)UA 5 Ufi+(5), 
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where 

L 

A* = ^(0) n U (£1 ° ^-r 1 ((-M)- 
3=1 

Here the sets Q±(5) are disjoint, compact, with finitely many connected components, 
while in general they are not invariant under the ifp-flow. 

Recall that F stands for the common value of the average of q over the tori Aj, 
1 < j < L. We introduce the following global assumption: 

There exist Ni, N 2 G N\{0} and a sequence 5 = 5j — > such that (1-24) 

inf ((Req) s -N 1 K -ReF)> 5 N \ sup ((Req) s -N 1 K - ReF) < -5 N2 . 

Here < K e is as in f)1.17[) . and we adopt the convention that when — 1, 
±5 N2 in the right hand side of (jl.24|) should be replaced by ±5/Ci for some C\ > 0. 

Theorem 1.1 Let aij and a 2 j be the fundamental cycles in Aj, 1 < j < L, defined 
by 

«ooj(afej) = {x e T 2 ; x k = 0}, k = 1, 2. 

VKe write then Sj = (Sij,S2j) and kj = (k(aij), k(a2j)) for the values of the 
actions and the Maslov indices of the cycles, respectively. Let us make the global 
dynamical assumption (jl.24jl . and assume that the differentials of the functions Pooj 
and Re (qj) , defined in (jl.22j) and (ll.23jl are linearly independent when C, = 0, 1 < 
j < L. Assume furthermore that e = 0(h 5 ), 5 > 0, satisfies e > h K , for some K 
fixed but arbitrarily large. Let C > be sufficiently large. Then the eigenvalues of 
P e in the rectangle 

h 5 eh 5 
\Rez\ < — , |Imz-eReF|<— (1.25) 

are given by 

P^( h (k-^-^ q h\+0(h^ kez\ \<j<L. 

Here Pj°°\^, e; h) is smooth in £ G neigh(0, R 2 ) and e G neigh(0, R) ; real-valued for 
e = ; and has an asymptotic expansion in the space of such functions, 

oo 

P}~\Z,c,h)~^h l p$\Z,e), l<j<L. 

1=0 

We have 

,(»)/ 



P^(e,e)=Poc,i(0+^)(e) + O(e 2 ). 
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Our next result treats the case when the strength of the perturbation e is suf- 
ficiently small but independent of h. In this case, we obtain a complete spectral 
result in a fixed /i-independent domain. 

Theorem 1.2 Let us continue to write Sj and kj for the actions and Maslov indices 
of the fundamental cycles in Aj, 1 < j < L. Assume that h 1 ^ 3 ^ 5 < e < €q <C 1, 
for some 5 > 0. As in Theorem 1.1, we make the assumption (I1.24J1 and assume 
that the differentials o/Poo,j(0 an d R e (<7i)(£) are linearly independent for I; = 0, 
1 < j < L. Let C > be large enough. Then the eigenvalues of P e in 



\Rez\ < 



A/N 



Im z 



ReF 



< 



A/N 



;i.26) 



are given by 
with 



n=0 



hik-h 



^,e) k e Z-'. I <./<!. 



holomorphic for £ = 0(e l l N ), and 



n 



0,1,2, 



1 < j < L, 



:-l 



Here pj is real on the real domain and the differentials of Pj{£,) and Regj(£,e) are 
linearly independent when £ = e = 0, 1 < j < L. We have 



The parameters N and N/ N can be taken arbitrarily large. 



a = aj. 



Remark. As we shall see in section 2, the assumption (jl.24j) implies the existence 
of a suitable weight function, which allows us to microlocalize the spectral problem 
for P t in a rectangle of the form (jl.25|) or (jl.26j) to a small neighborhood of the 
union of the tori A^, 1 < j < L — see Proposition 2.3. Indeed, it is the existence 
of the weight function that allows us to carry out the complete spectral analysis 
in such rectangles. The purpose of the condition ()1.24|) is to provide an explicit 
criterion of a purely dynamical nature, which suffices for the construction of the 
global weight. As will also be seen in section 2, in the case when the if p -flow is 
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completely integrable, the condition ()1.24|) can be replaced by a slightly different, 
although equivalent one, having the advantage of being easier to verify in practice — 
see (|2.53|) below, Proposition 2.5, and the discussion in section 7. In particular, in 
section 7, we show that in the completely integrable case, the set of values F in 
(jl.21j) to which Theorem 1.1 and Theorem 1.2 apply, covers the entire region ()1.19|) . 
apart from a subset of a suitably small measure. Such a result is stated in Theorem 
7.6, where, using the isoenergetic KAM theorem, we also extend it to the perturbed 
situation, when the completely integrable symbol p is replaced by p + O(X), with 
A > small enough. 

The construction of quasimodes associated to invariant Diophantine Lagrangian 
tori and to Cantor families of such tori has a long tradition in the self adjoint case, 
and we refer to [IHj, [Ej, [2Z|, and |3Z| for the results in this direction. Especially 
relevant in the present two-dimensional context is the work of Shnirelman jHEl m the 
selfadjoint case, which contains the idea of using invariant Lagrangian tori to split up 
the real energy surface into different invariant regions. See also [TI]. In jSH], the main 
focus is on constructing quasimodes associated with the gaps between the invariant 
tori in a perturbative situation. In our non-self adjoint case, the idea of using the 
invariant tori as separatrices becomes more efficient than in the standard selfadjoint 
setting, and in a future work we plan to use it to study the global distribution of 
eigenvalues inside the entire band corresponding to different invariant tori. 

In the present work, we essentially exploit only the invariant tori ()1.13|) . corre- 
sponding to a given value of the average of the leading perturbation, and obtain 
complete spectral results in the spirit of our previous works jH], [22], and [21]. In 
particular, the work [21] introduced and exploited a dynamical condition somewhat 
reminiscent of (jl.24j) . Technically speaking, many important differences occur, how- 
ever, due to the more general behavior of the classical flow away from the union of 
the invariant tori (|1.13|) . 

The plan of the paper is as follows. 

In section 2, we use the global dynamical condition (jl.24|) to construct a globally 
defined compactly supported weight function, which allows us to microlocalize the 
spectral problem for P e to a small neighborhood of the union of the Ay's, 1 < j < 
L. We also give a discussion of a modified version of ()1.24|) . when the if p -flow is 
completely integrable. 

In section 3, we carry out a construction of a quantum Birkhoff normal form for 
P e , valid to an arbitrarily high order along the invariant torus, and in (e, h). 

In section 4, we solve an appropriate Hamilton- Jacobi equation, to be used in 
the spectral analysis of P t in the case when e is sufficiently small but independent 
of h. The ideas used here are similar to [2E] and 0T] . 
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In section 5, we justify the eigenvalue computation based on the Birkhoff normal 
form construction from section 3 and prove Theorem 1.1, by solving a suitable global 
Grushin problem. 

In section 6, we use the results of section 4 to modify the Grushin analysis from 
the previous section, in order to establish Theorem 1.2. 

In section 7, we first study the completely integrable case and verify the global 
dynamical condition (jl.24j) . or rather (j2.53|) below, under some general assump- 
tions. This is then applied to the case of convex analytic surfaces of revolution, and 
complex perturbations, close to rotationally symmetric ones. We then exploit the 
isoenergetic KAM theorem and study the dynamical condition in the case when the 
unperturbed symbol is close to a completely integrable one. The results of this sec- 
tion are summarized in Theorem 7.6, which together with Theorem 1.1 and Theorem 
1.2, can be considered as the final result of the present work. 

In section 8, we give an application of Theorem 1.2 to the barrier top resonances 
for the semiclassical Schrodinger operator and obtain an extension of the result 
of j2E] to an /^-independent region in the complex plane. 

Acknowledgment. Part of this project was conducted when the first author vi- 
sited Ecole Polytechnique in June of 2004. He is happy to thank its Centre de 
Mathematiques for a generous hospitality and excellent working conditions. The 
partial support of his research by the National Science Foundation under grant 
DMS-0304970 is also gratefully acknowledged. All three authors have benefited from 
the hospitality of the MSRI, Berkeley, in Spring 2003, during a special semester on 
semiclassical analysis. 

2 Construction of the global weight 

In and [21] , the basic weight function was coming from an averaging procedure 
along the if p -flow, and was given by Gt, defined in (I1.18|) . with T being a common 
period for the flow. This weight was employed in a full neighborhood of p _1 (0) H 
T*M. As we shall see in section 3, in the present case, the Birkhoff normal form 
construction will in general be valid only to an infinite order along an invariant 
torus, which will force us to work in small /i-dependent neighborhoods of the tori. 
We shall therefore find it convenient to reduce q to the flow average {q)r,K by means 
of Gt when away from a small but fixed neighborhood of the union of Aj, 1 < j < L, 
while near each Aj, we shall use a somewhat different sort of average, which is close 
to (q)r,K, when T is large. 

We keep all the assumptions of the introduction, and consider an operator P t 
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with the leading symbol 

p + ieq + 0{e 2 ), (2.1) 

microlocally in a neighborhood of p _1 (0) n T*M. 

Let us say that a smooth multi-valued function defined on or near a torus is grad- 
periodic if its gradient is single-valued. In what follows, we shall work microlocally 
near a fixed invariant torus, say Ai. From the assumption ()1.13j) . we recall that 
there exist analytic and grad-periodic coordinates on A 1; x\ and x%, which induce 
an identification between Ai and T 2 . Applying Weinstein's tubular neighborhood 
theorem, see j^], or simply following the argument, described in the beginning of 
section 1 of [2H], we see that we can extend X\ and x 2 to a tubular neighborhood 
of Ai and complete them with analytic functions £1 and ^> defined near Ai and 
vanishing on this set, so that we get a real analytic canonical transformation, 

K\ : neigh(A x , T*M) -> neigh(£ = 0, T*T 2 ), (2.2) 

which maps Ai to the zero section in T*T 2 , and such that when expressed in terms 
of the coordinates x and £, the unperturbed leading symbol p of P e=0 becomes 

p(x,0 = a-£ + O(e), (2.3) 

with a satisfying ()1.15|) . It will now be convenient to perform an additional real 
canonical transformation on T*T 2 , in order to make p independent of x to a high 
order in £. To this end, we observe that a straightforward Birkhoff normal form 
construction, very similar to the one described in detail in section 3, shows that there 
exists a sequence of real- valued functions G\,G 2 ■ ■ ■ , with Gj being homogeneous of 
degree j + 1 in £ and depending analytically on x, such that if 

G ~ G\ + G 2 + • • • , 
then at the level of formal Taylor expansions in £, we have 

poexp(H G )=p oo (0 + O(D, 

where Poo(C) = a ' £ + C(£ 2 ) depends on £ only. See also (jl.22|) . Since we work 
in the analytic category and do not wish to consider convergence questions for the 
Birkhoff normal forms and associated canonical transformations (in this connection, 
see [30J), we truncate the series at some fixed but arbitrarily large order N, and 
write 

poexp(H G(N) )=p N (0 + O(Z N+1 ), 
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with Cr(jv) = Gi + G 2 + ■ ■ ■ + Gjf-i, and Pn(0 = a • £ + C(£ 2 ) independent of x. 
Notice that since Grm(£,) = 0(^ 2 ), the analytic canonical transformation 

cap (ff 0(W) ) (2.4) 

maps the zero section £ = to itself and it also preserves the action integrals along 
closed loops. The composition of the transforms, with K\ defined in ()2.2|) . 

maps symplectically a neighborhood of Ai in T*M onto a neighborhood of £ = in 
T*T 2 , in such a way that Ai is mapped onto the zero section £ = 0. Implementing 

KN l = K ( N ) o Kl (2.5) 

by means of a microlocally unitary Fourier integral operator with a real phase, U, 
and conjugating P e by means of U, we obtain a new operator, still denoted by P e , 
which is microlocally defined near the zero section in T*T 2 , and has the leading 
symbol 

p (x, £, e) = p(x, + teq{x, f) + 0(e 2 ), (2.6) 

with 

p(x,0=Piv(0 + O(^ +1 ), PN(0 = a-Z + O(e)- (2.7) 

Here the operator P e =o is selfadjoint. Moreover, the complete symbol of P e is a 
holomorphic function in a fixed complex neighborhood of £ = 0, and it depends 
smoothly on e G neigh(0, R). On the operator level, P e acts in the space of microlo- 
cally defined Floquet periodic functions on T 2 , Lg(T 2 ) C L 2 oc (R 2 ), elements u of 
which satisfy 

u(x -v) = e ie - u u(x), 6 = ^- + —, v G 2vrZ 2 . (2.8) 

2%h 4 

Here S = (S±, S2) are the classical actions, 

Sj = / £dx, j = 1,2, 

Jotj 

and otj form a system of fundamental cycles in Ai, defined by 

Kite) = ft, i = 1,2, Pj = {xe T 2 ; = 0}. 

The tuple ko = (fcote), ^ote)) stands for the Maslov indices of the cycles ay, 
J = 1,2. 
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In what follows we shall work with the operator P e , microlocally defined near £ = 
in T*T 2 , which has the principal symbol ([2.60 . (|2.7|) . For future reference, let us 
also remark that an application of the implicit function theorem shows that we may 
assume that the level set pn = is of the form 

6 = /(&), (2-9) 
for some analytic function / with /(0) = 0, /'(0) ^ 0. 

We shall now discuss the problem of solving the equation 

H p G = q-r, (2.10) 

considered near £ = in T*T 2 , for a suitable remainder r. The torus average of the 
left hand side of ()2.1U|) is 0(£ N ), and we shall try to make r independent of x at 
least to that order in £. Also, the function G should be analytic near £ = 0. When 
solving (|2.10|) . we first consider the problem of solving 

H PN G = q-r. (2.11) 

We analyze this problem on the level of formal Taylor series in £, and look for G in 
terms of a formal expansion, 

oo 

g ~ y^G^, 

fc=0 

where is homogeneous of degree k in £. Taking also a (finite) Taylor expansion 
of with respect to £, we easily see that 

oo 

H PN G ~ 2J /«> 

n=0 

where /„ is homogeneous of degree n in £ and is given by 

fn = {PV,Z, Gfc}, 

fc+i-l=n 

with j?jv,z homogeneous of degree / > 1 in £. Introducing also a Taylor expansion 
of q with respect to £, q ~ go + 9i + • • • , an d using the fact that the operator 
a ■ d x is globally analytic hypoelliptic, we determine successively analytic functions 
Gq,Gx, . . ., with Gk homogeneous of degree k in £, such that 

q n - fn = (q n ), n = 0, 1, . . .. 
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Here, as in (|1.23p . we write 



(/>(0 = T^p / f&Odx, 

for the torus average of a smooth function / defined near £ = in T*T 2 . Taking 
the finite sum Gq + G\ + . . . + Gn, we get an analytic function G such that 

H PN G = q-r, 

where 

To solve ([2.10)1 . we just use that 

H p G = H PN G + Ofc N ). 

We conclude therefore that for any fixed N e N as in (|2.7|) . we can find an analytic 
function G defined in a fixed iV-independent neighborhood of £ = and solving 
(|2.10|) . with the remainder r satisfying 

r(x,0 = (q)(0 + O(Z N )- (2.12) 

Composing G with the inverse of «jv,i, defined in (|2.5jh we get an analytic function 

Gi^Go^, (2.13) 

defined in a small but fixed neighborhood of Ai C T*M. 

Coming back to the torus model, let us recall that we are working in real symplec- 
tic coordinates (x, £) for which ()2.7|) holds. We shall consider the behavior of the H p - 
flow near £ = for large, but finite times. With (x(t),£(t)) = exp(tH p )(x(0),£(0)) 
we have by the Hamilton equations, 

m = o(\m N+1 ), 

and a standard " continuous induction" argument shows that for \t\ < 0^(1) |£(0)| N . 
we have 

m = o(\m\), 

provided that |£(0)| is small enough. We get, for these times, 

m = m+to(\m\ N+1 )- (2.14) 
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Similarly, we have 

i(t) = ^r+mn (2.i5) 

and with (ar o (t),£(0)) = exp (tH PN ){x{0), £(0)), we get for |t| < 0^(1)^(0)1"^, and 
writing (x,£) rather than (x(0), £(0)), 

e(t)=e + ^(Kr +1 ), x{t) = x {t)+tO{\t\ N )- (2-16) 

For future reference, we shall also investigate the behavior of (q)T,K,p near £ = 0, as 
T — > oo. In doing so, we notice that ()2.16|) gives 

(«>Wa?,0 = <?>T^ w (*,0 + O(T|eH, r<O w (l)|e| _JV . (2.17) 

where we may also remark that the 0-term in ()2.17j) depends only on the Lipschitz 
norm of q. When analyzing (q)T,K, PN , we write 

and expanding g(-,0 in a Fourier series, we get 

(q}T,K, PN &0 = (l)(0+ E e <ta g(fc,0^(Ti/ w (0-*). (2-18) 

o^fcez 2 

Here <f(fc,£) are the Fourier coefficients of g(a;, £) and = / e" lt ^K{t) dt is the 

Fourier transform of f( . Now the Diophantine condition ()1.15|) gives 

IpUO ■ *i > — ^ - d iei 1*1 , mo, c 1 >o, 

G |fc| 

and therefore 

\Pn(0 ■ k\ > 1 -^ r , if C x \i\ \k\ < 1— jr , fc ^ 0. 

^ 1 20,1*1° 2C l/fcl^ 

Let < x G C 3O ((— 1, 1)) be equal to 1 on (—1/2, 1/2), and let us decompose the 
sum in (I2.18J) as follows: 



o^fcez 2 

+ E (2GA |£| |Af° +1 ) e fa '*g(fc, 0^(^(0 • *) = I + 11, 



o^kez 2 
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with the natural definitions of I and II. We get for each M G N, 
X (2CcA |f | \k\ No+1 ) 

\i\ < 0m(i)E , t ; m ttm 7 1^01 < Cm(1)T- m ^|A;| MJVo 

and therefore, 

/ = C(T-°°). 

When treating II, we only use that K is bounded and because of the presence of 
1 — x, we have \k\ > |f \~ 1 /( No+1 * in the sum, so that the contribution coming 

from II is 0(\£ |°°). We conclude that 

(q) T ,K, PN &0 = (q}(0 + O(C + T-°°)- (2-19) 
The discussion above is summarized in the following proposition. 

Proposition 2.1 Let (x, f) be real symplectic coordinates so that ()2.7|) is true. Here 
N is fixed but arbitrarily large. We then have for f small enough, 

(q)T,K, P (^0 = (q)(0 + O(r + T-°°)+O(T\Z\ N ), T < O n {1) \£\~ N . (2.20) 

For future reference, we shall pause here to recall some general estimates for 
convolutions of the form Kt * g, where g G L°°(R). The following discussion is 
motivated by the fact that the flow average of q, 

(q)T, K (p)= J K T (-t)q(exp {tH p ){p))dt 

is a convolution in the time variable along the if p -trajectory passing through p. The 
starting point is that 

K = K e * K + r e , ||r e || L i = 0(e), (2.21) 
for < e < 1. It follows that for < S < T, 

K T = K s * K T + rs T , (2.22) 



where rs T (t) = ^rs(-f), so that 



rs^ T \\ L i = || rs \\ L i = £>(-). 



S 
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rs T *g\\ L «,<0[^)\\g\\ L oo. (2.24) 



If g G L°°(R) is real- valued, we have 

Kt * g = Kt * Ks * g + rs T * g, (2.23) 

with 

1 1 r .<? „ * n 1 1 r no < D I 
In particular, when || g || £ oo = 0(1), we get 

inf (K T * g) > inf (K s * g) - O ) , (2.25) 

and 

sup(K T * g) < sup(K s * g ) + o(~). (2.26) 

R R \ -L 



Here we can be a little more precise about where to take the sup and the inf. If 
I C R is an interval, then 

^(K T *g)>^M pK (Ks*g)-0^), (2.27) 



sup{K T *g)< sup (K s *g) + 0[^). (2.28) 

I I-TsuppK \ 1 



We also notice that if K has the same properties as K, then (|2.21|) can be generalized 
to 

K-K*K e = 0(e) inL 1 , 
leading to the possibility of replacing K s by K s in (l2~22l . (12~23D . and (12~231) - (12~2%|) . 

Let us return now to the operator P e with the leading symbol ()2.6|) . ()2.7|) . We 
shall assume from now on that 

rfpiv(O) = a and <iRe (g)(0) are linearly independent. (2.29) 

Notice that this condition is independent of N. Then, possibly after changing the 
sign of the ^-coordinate and using ([2.9)1 . we get 



_d_ 



(R*<7>(6,/(£i))>0. (2.30) 

6=o 
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Lemma 2.2 Assume that ()2.29j) holds true, and let N in ()2.7|) be sufficiently large. 
There exists a constant Cq > such that if (ll.24j) holds for some 1 < Ni,N 2 , and 
Ni G N\{0} ; then there exists 5q > such that 

inf ((Req) s _ Rl K - ReF) > — , < 6 < 6 , (2.31) 

and 

sup f<Reg)._^ „ - ReF) < - — , < 5 < 5 . (2.32) 

H_(<5) v ' /O 



Proof: It suffices to prove only (J2.31|) . and we may also assume that F = 0. More- 
over, as will be clear, the following argument will not depend on the choice of N in 
(|2.7p . provided that it is large enough, and to simplify the discussion notationally, 
we shall take N = oo. We shall work near one of the tori Aj. From Proposition 2.1 
we then know that for every fixed N\ G N\{0}, we have 

(g>^x^(x > o = <ff>(o+o(n. iei<^, 

if < 5 < 5(Ni) > 0. Take 5 as in and let 5 G [5 Ml ,S M °], where 1 < M < M x 

will be fixed later. Then 

inf (Rfig)j_jf > min inf (Reg)^ inf (Reg) , (2.33) 

n+(<5) ' yn+(5)\n+(25) ' ^+(2<5) ^ J 

and we also get 

Jnf (Iteg)-.^ >_ (2.34) 

n + (5)\n+(2<5) ' o 

for some Cq > 0, when 5 is small enough. We may assume that Co > 2. 

If we choose M with M A^x > N u then 5 _JVl > 8' Nl , and using the general 
estimate ()2.27|) in the form 

(q) T ,M > t m^(q) s ^w(tH p )(p)) - O 

and the fact that in the region |£| < 5, the if p -trajectories are approximately confined 
to the tori £ = Const, for times bounded by any fixed inverse power of |£|, we get 

inf (Re g)^ K > inf (Reg) - O(r^) > ^ - O(r^). (2.35) 

S2_l_(2oJ ' S Z_|_ (oj 



18 



Here we should replace 5 N2 by 8/C\ if N 2 = 1. Now assume that M > N 2 , M Ni — 
Ni > N 2 , so that 5 Na > 8, 8~ Nl 8^ < 8 N \ Then fl2~55j) implies that 



Combining this with ()2.33|) and (j2.H4jl . we get 

inf (Reg)^ > — , 

n+(S) ' 

for 5 > sufficiently small and for 8 G [5 A/l , 5 M °], provided that 

N 1 + N 2 



Mi > M > max iV 2 



iVi 



and 5 > is small enough depending on the exponents and such that the estimate 
flOU holds true. 

We conclude that we have the estimate in (jl.24j) with N% = Ni and N 2 = 1 for 
S replaced by 8 G [8 Ml ,8 °]. A special case of this is when we start with 8 with 
Ni = Ni, N 2 = 1, and then conclude that (|1.24|) holds with 8 replaced by 8 in 
[8 Ml ,8 Mo ], when M 1 > M > + 1)/A\. Choose M x = M 2 . Then the argument 
can be iterated and we get (jl.24j) with A^i = Ni, N 2 = 1, for 8 replaced by any 8 in 
[8 M °, 8 Mo ) U [8 M o, 8 M °] U . . . = (0, 8 Mo ), and the lemma follows. □ 

Lemma 2.2 shows that (jl.24)l is independent of the choice of Ni and N 2 . The 
proof also shows that it is independent of K G Cq°(R). We also notice that the 
proof shows that the validity of (jl.24|) is stable under small perturbations of p that 
preserve the invariant tori Aj and the Diophantine condition (jl.l5j) . 

We shall now construct a global weight function G. When doing so, let us recall 
from (J2.13|) the analytic functions Gj = G o Kjf-, defined in small neighborhoods of 
A-,-, 1 < j < L. Here k n ^ is defined exactly as in (|2.5|) . When constructing 
the global weight G, we shall glue together the analytic functions Gt, defined in 
(J1.18)) and Gj. In doing so, it will be sufficient to work in a fixed sufficiently small 
neighborhood of f = in T*T 2 . Let < x < h X e C^(R 2 ), be supported in 
|£| < 2, with x = 1 on |^| < 1. With < fi <^ 1 to be determined, we put 



G= t _y ) G T + Yv G . (2 .3 6 ) 
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In the following calculation, in order to ease the notation, we shall take L = 1, in 
which case we may omit kjv,i from the notation and write 

G= (I-X^Gt+X^G, (2.37) 

where = x(C/ 7-0 ■ We shall compute H p G, and in doing so we notice first that 

H pX , = Oft"), 

uniformly in fi > 0. We get 

H P G = (1 - (? " (q)T,K, P ) +XAQ- (?» + C(e JV )XM + # p Xm (g - G T ) , 
which we rewrite as 

H P G = q - (q) - (1 - X,)((q)T,K, P - (g» + + - <?t). (2.38) 

Now x M = O(l) and fTTTHJ) gives G - G T = 0(T), so that 

H P G = q - (q) - (1 - XmX^t,^ - (g» + 0(T£ N ), T > 1. 
Proposition 2.1 gives then that uniformly in /x > we have 

= 9 - (q) - (1 - X/1 ) (O(T^) + 0(r + T- 00 )) + 0(T£ N ), T < O n (1) \£\- N . 

(2.39) 

We let first N be sufficiently large depending on Ni, and take then 5 > small 
enough but fixed. Take then /j ~ 5, and put 



T = li 



-Ni 



In the region |£| < /i, we have G = G, while G = Gt for |£| > 2/x. In the intermediate 
region where < Xn < 1> we have |£| ~ /i, and it follows that in this region 

q -H p G = (q)(0 + O(Z N - Nl )- 
The discussion above can be summarized in the following proposition. 

Proposition 2.3 Assume ()2.29|) for each 1 < j < L ; and (jl.24|) /or some iVj. and 
iV 2 G N\{0}. When 1 < j < L ; letGj be an analytic solution near A j of the equation 
(|2.1(J|) . mi/i r satisfying (|2.12|) /or some iV Zarge enough. Let Wj be a sufficiently 
small neighborhood of Aj, Wj = (£i o /Civj) -1 ((— <5 , 5q)) ; 5o > 0. Then there exists 
a smooth compactly supported function G on all ofT*M, such that G = Gj in Wj, 
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1 < j < L, and if W = lif =1 Wj, 5{W) > is sufficiently small, and Q±(5(W)) 
stand for the components of Q±(Sq) in p~ l ((— S(W), 5(W))) \W , we have 

Re( 9 -/f p G-F))>^, 

Re( 9 - Hp G-F))<^. 
Here F is the average of q over Aj, 1 < j < L, and C(W) > 0. 

We shall now discuss a modified version of the global assumption (jl.24j) . available 
if we assume that the if p -flow is completely integrable. In particular, the results 
of the following discussion will be applied in section 7, where, as an example, we 
consider p corresponding to the geodesic flow on an analytic surface of revolution. 

Assuming the complete integrability for the flow, we know that there exists an 
analytic real- valued function /, which Poisson commutes with p, and such that each 
Aj, 1 < j < L, is a level set for the associated mapping, 

(p, /) : T*M -> R 2 , (2.40) 

corresponding to a regular value. Then each torus Aj is embedded in a Lagrangian 
foliation of if p -invariant tori, given by 

A E , F :p = E, f = F, 

for (E,F) G neigh((0, Fj), R 2 ), where {0,Fj) is a regular value of ()2.40|) . We may 
then introduce the action-angle coordinates, given by a real analytic canonical trans- 
formation, 

Kj : neigh(A i; T*M) -> neigh(£ = 0, T*T 2 ), 1 < j < L, 

such that Aj is mapped to the zero section in T*T 2 , and such that when expressed 
in terms of the coordinates x and £, the unperturbed leading symbol p becomes a 
function of £ only, with 

p(Z)=a-£ + 0(e), a = a v (2.41) 

As in (12. 9j) . we may assume that the energy surface p _1 (0) takes the form £2 = 
/.,(£i), where fj is analytic with /j(0) = 0. The tori A 0) f, for F G neigh(Fj,R) 



and 



inf 

ci+(6(w)) 



sup 

Q-(8(W)) 
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then take the form ^ = /i, £ 2 = fj(fj), for < b, with < 6 <C 1. Consider a 
flow-invariant neighborhood of the union of the Aj of the form 

L 

U U(eioK i r 1 (A*)n P - i (o)= u a,. 

\n\<bj=l H<6 

Here is a disjoint union of L flow-invariant tori. We decompose the real energy 
surface as follows, 

p-\0) = A_ b U [J A,UA 6 , 

where A± b are two disjoint compact flow-invariant domains with at most finitely 
many connected components. This gives a decomposition of the energy surface 

p-\0) = |J A„, (2.42) 

where M — [—b, b] and each A M is a compact if^-invariant set, with finitely many 
connected components A^j. Assume also the continuity property: 

For any /i G M, e > there exists 5 > 0, such that for (2.43) 
dist(//,/i ) < 8 we have A M C A M0i6 := {p G p _1 (0); dist (p, A M0 ) < e.} 



As in (jl.l7j) . we put 

(q)t,M 



q(exp(tH p )(p))K T (-t) dt. 



(2.44) 



Let Qr{^ij) and Qr(fj) be the ranges of Re (q)r,K restricted to A M j and A M , respec- 
tively. Then Qt(^) is a finite union of the closed intervals Qt(h, j), and according 
to (ET2TI) and (l2~2KJ) . we have 

T'T 



Q T (p)cg 5 (/i) + o(i)lkl 



L°° 



(2.45) 



uniformly for < S < T, \l G M. Let Qoo(p, j) be the non-empty intersection of 
all the Qt{^i j) for T > 1, and define Qoo(p) to be the union of all the Qoo(p, j), so 
that 



Qoo(^) = [j 



lim inf Re {q)T,K, lim sup Re {q)r,K 



T— >oo A 



(2.46) 



We may also remark that according to (|2.45|) . Qoo(n) C Qt (/•*), for all T > 1. It is 
also easy to see that <3oo(/-t) does not depend on the choice of K. Put 

Q T = {(//, E); ii G M, £ G Q T (/x)}, < T < oo. (2.47) 
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Lemma 2.4 We have 

1. Qt is closed. 

2. For every neighborhood U of there exists a T G (0, oo) such that Q T C U 
forT>T . 

Proof: We prove (1) first. It is clear that Qt is closed for each finite T. Let 
(fjLj, Ej) G Qoo be a convergent sequence so that jij — * /to G M, Ej — > i?o G R. Let 
e > 0. Then there exists T G (0, oo) such that Qt (^o) C Qoo(a*o) + [ — e ) e ]- We fix 
such a number. If 5 > is small enough, we have 

QooGu) C Q To (/i) C Qt (//o) + [~e,e] C Qoofao) + [-2e,2e], (2.48) 

for dist(/i, /io) < 5- In particular, _E G Qoo(^o) + [— 2e, 2e], and letting e — > 0, we 
get E'o G Qcc(n>o), proving the closedness of Qoo. 

We now establish (2). If T = kT , k > 1, we get from (jO^) and (fSHHD that for 
dist (/i,/i ) < 5, 



Qr(A*) C Q^o) + [-2e,2e] + 

and hence 



0(1) 0(1) 



k 1 k 



Qr(fJ-) C Qoo(Aio) + [-3e, 3e], 

for T large enough. This means that {(//, £?) G Qt; dist (/x, /to) < 6} is within a 
distance 3e + 5 from Q M . The statement (2) now follows by a covering argument. 
□ 

If |/t| < b, then each connected component A w -, 1 < j < L, of A M is diffeomorphic 
to T 2 in such a way that H p \\ . becomes a\d Xl + c^c^, aj G R, and we have 

fcez 2 

which, as T — > oo, converges uniformly in x (using now also the smoothness of q), 
to 

(?WaO= Eg(fc,/x)e^' fc . (2.50) 

a-fc=0 

Notice that Qooil^ij) is the range of this function. 

When |/t| < b and 1 < j < L, we shall now study the compact interval Qooifaj) 
in more detail. In doing so, we remark that if we have an equation of the form (|2.1(J|) 
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on U^^Apj, with G smooth and bounded, then Qoc(n,j) C Rer(A At ). In fact, with 
K = l[-i,o], we have for p G A MJ -, 

(r) T (p) = (q) T (p) ~ f (G(exp (W p )(p)) - G{p)) = (q) T (p) + O 

and it suffices to let T —>• oo, and use that Re (r) T (A M ) C Rer(A /J ). 

Recall now that we have seen that we can solve (|2.1(J|) with r satisfying (J2.12j) . 
and we conclude that 

Qoo(v,j) C Re (<&)(//, /(//)) + 0(ji N )[-l, 1], < b, q 3 = q. (2.51) 

As in the general case, from now on we introduce the assumption (J2.29|) for each 
1 < j < L. In view of ()2.51|) . this implies (possibly after shrinking b again) that for 
< b < b, 

inf |J Q M (/i) >ReF + ^, sup |J < Re F - A (2.52) 

A»S[6,6) /Lte(-6,-6] 

for some constant C > 0. The global assumption in the completely integrable case 
is then the following one: For any b G (0, b) there exists C(b) > such that 

inf QooifJ.) > ReF + — L-, sup Qoo(aO < ReF l —. (2.53) 

M 6[6,6] C(b) pe[-b,-b] C{b) 

Using (|2.51|) together with Proposition 2.1 and Lemma 2.4 it is easy to see that 
in the completely integrable case, the condition ()2.53|) is equivalent to () 1.24)1 . 

Proposition 2.5 Assume that the H p -flow is completely integrable, so that ()2.41)) 
holds true. Assume furthermore that (j2.29j) and ()2.53j) are valid. Then we have the 
same conclusion as in Proposition 2.3. 

In what follows, to ease the notation, we shall drop the tilde and write G instead 
of G, defined in Proposition 2.3. Associated with G G C^°(T*M), there is a globally 
defined IR-manifold 

A eG = {p + teH G (p); p G T*M} c T*M, (2.54) 

and when acting on the Hilbert space H(A e a), associated to A e c by means of the 
FBI-Bargmann transform 

Tu{x) = Ch- 3/2 / e V h V u{y) dy, (p(x, y) = -(a; - y) 2 , 
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the operator P e gets the leading symbol 



(p + ieq + 0{e 2 )) {p + ieH G p) = p{p) +ie{q- H p G) (p) + 0{e 2 ). (2.55) 

Here we are tacitly assuming that M = R 2 so that T*M is a linear space. In the 
manifold case, using that G is analytic near Aj, 1 < j < L, we define A eG = 
exp (ieHc)(T*M) in a complex neighborhood of Aj, and elsewhere we let G G 
C^(T*M) stand for an almost holomorphic extension of G, and A e c is then de- 
fined as 

A eG = exp(eff^2)(r*M). 

Here the Hamilton vector field of Re G is computed^with respect to the imaginary 
part Im<7 of the complex symplectic form a on T*M. 

It follows from Proposition 2.3 that when m e A eG is away from a small but 
fixed neighborhood of the union of the Lagrangian tori Aj, such that each Aj is the 
image of Aj in A eG under the map T*M3p*->m = p + ieHa(p), it is true that 

llm P e (m) - eReF | > 



o(iy 

provided that |ReP e (m)| < 1/C for a sufficiently large C > 0. Moreover, microlo- 
cally near p _1 (0), the operator P e acting on H(A e c) is unitarily equivalent to an 
operator acting on L 2 (M), which has the leading symbol (|2.55|) . 

Summing up the discussion of this section, we have achieved that microlocally 
near each Aj C A e c, I < j < L, the operator 

P e : H(A eG ) -> H(A eG ) 

is unitarily equivalent to an operator P e , acting on L^(T 2 ) and defined microlocally 
near £ = in T*T 2 . Here P e is such that 

oo 

P e ~ y *r / h''p v (x,Z,e), (2.56) 

with p v holomorphic in a fixed complex neighborhood of £ = and 

Po = Pn(0 + ie(q)(t) + 0{e 2 ) + 0(^ +1 ) + eO(^), (2.57) 

with Pn{£) = cl • £ + 0{^ 2 ), and where the term (9(£ 7V+1 ) is real on the real domain. 
Here is fixed but arbitrarily large. In what follows, we shall drop the tildes and 
write P t and p v , v > 0, instead of P t and p v , v > 0, respectively. 
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3 The normal form construction 



We recall that we have reduced the analysis to the operator P e with a complete 
symbol ()2.56|) and a leading symbol (|2.57|) . and recall also the assumption (|2.29p . 
Our goal in this section is to construct a quantum Birkhoff normal form for P e — see 
also HOI, [S2], and jaZJ- We may write 

Po(x, £, e) = po^x, £, e) + p , 2 (x, £, e) + . . . + p , N (x, £, e) + 0((e, (3.1) 

where p 0) j is homogeneous of degree j in (£, e), so that in particular 

Po,i(x,£,e) = a-£ + ie(<?}(0) 

is independent of x. 

We shall now remove the x-dependence also in the terms P qj, 2 < j < N. In 
doing so, we consider the formal power series 

oo 

G = J2G J (x,£,e), (3.2) 

3=1 

where Gj = 0((£, e)- 7 " 1-1 ) depends analytically on x and is a homogeneous polynomial 
of degree j + 1 in (£, e). In the sense of formal Taylor expansions, we then have 

oo oo oo 

Po o exp (H G ) = po + & H gPo = + fcT^Po,! 

fe=l ' fe=l Z=l 

oo oo oo oo ^ oo 

fc=l i=l ji=l j k =l ' n=l 

Here q n is a homogeneous polynomial of degree n in (£, e). We also notice that in 
the sum above Hc m Po,i is homogeneous of degree m + I in (£, e). It follows that 
q x = Po,i is independent of x, g 2 = Po,2 + H Gl p 0:1 , and g n+ i = p ,n+i + H Gn p 0:1 + q n+u 
n > 2, where depends only on Gi, . . . G n -\. We are therefore able to determine 
Gi, Gi . . . Gn-i successively by solving the cohomo logical equations 

H P0,l G n = a ' d xG n = p ,n+l + 9n+l ~ (Po,n+l + 9n+l), 1 < Ti < iV — 1. (3.3) 

In this way we achieve that all the q n , n < N, are independent of x. 
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Remark. It follows from the construction that the function G\ is independent of £ 
and is of the form Gi(x,£,e) = e 2 h(x) for some analytic function h. Furthermore, 
it is clear that 

G j = O j (e 2 ), 1<3<N-1. 

Remark. In the construction above, we could also have used a slightly different 
argument and looked for G as a formal power series in e alone, 

G~G 1 (x,0 + eG 2 (z,0 + ... 

Then composing the symbol po with the canonical transformation exp (€ 2 Hg) and 
repeating the previous arguments, we would have been able to determine each Gj 
modulo 0(£ N+1 ), and the final result would have been the same. 

Summarizing the discussion so far, we get the following result. 

Proposition 3.1 Let p (x,£,e) = p N {£) + ie(q)(£) + £>(e 2 ) + 0(£ N+1 ) + eO(£ N ) 
be an analytic function defined near £ = in T*T 2 , depending smoothly on e G 
neigh(0, R). Here N e N is fixed but can be taken arbitrarily large. Assume that 

p N {t) = a-i + 0{e), 

where a satisfies (jl.l5|) . Then we can find analytic functions Gi,G , 2,... ,Gn-i, 
with Gj(x,£,e) being a homogeneous polynomial of degree j + 1 in (£, e), such that 
Gj = Oj(e 2 ), 1 < j ' < N — 1, G\ is independent of £, and such that if 

GW = G 1 + G 2 + ... + G N - 1 , (3.4) 

then 

p o exp (H G(N ) ) = p {N) (f , e) + r N+1 (x,£,e). 

Here 

p (N) (^e)=a-Z + ze(q}(0) + O((Z,e) 2 ) 

is independent of x and rN + i(x,£,e) = 0((£, e) N+1 ). Writing p^(£,e) = p(£) + 
ieq(£,e), where p is real, we have 

d,£p(0) Re d^q(0, 0) are linearly independent. 

In section 5, we shall see that we can quantize the holomorphic canonical trans- 
formation 

K := exp (H G (n)) (3.5) 
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by means of an analytic elliptic Fourier integral operator in the complex domain. In 
this section we shall proceed somewhat formally, and carrying out the corresponding 
conjugation of P e , we may assume from now on that we are given an /i-pseudodiffe- 
rential operator, still denoted by P e , defined micro locally near £ = in T*T 2 , whose 
full symbol has a complete asymptotic expansion, 

P e (x, £, e; h) = p + hp 1 + h 2 p 2 + ..., (3.6) 

with all Pj = pj(x,C,,e) holomorphic in a fixed complex neighborhood of £ = 0, 
depending smoothly on e G neigh(0, R), and such that 

p (x^,e)=p^^,e) + O(^,e) N+1 ), ^(e ) e)=a.(+ l e( 9 >(O) + 0(((, e ) 2 ). (3.7) 

Our goal now is to make the lower order terms Pj, j > 1, in ()3.6|) independent of x, 
to a high order in £ and e. This will be achieved by means of the usual conjugation of 
P t by an elliptic pseudodifferential operator of the form exp (Q(x, hD x , e; h)), where 
Q is of order in h. Write Q(x, hD x , e; h) ~ g + hq\ + . . . . Then on the operator 
level, we have 

oo oo 1 

e Q Pee -Q =P £ + J2J2yh l {^Qfpi 
k=l 1=0 

= ^ + E E E • • • E 7>^ + - + - +fc iadfoj ... ladfej p h 

k=l l=0 j 1= j k =0 ' V / V / 

and on the symbol level we get 

oo 

e Qp t6 -Q „ h n s n . (3.8) 

n=0 

Here s = p , Si = p\ + iH po q , and for n > 1, the term s n+ i has the form = 
Pn+i + iH po q n + s n+ i, where depends only on q , . . . q n -i- We wish to make the 
terms s n , n < N, independent of x, modulo 0((£, e) Taking a Taylor expansion 
of go in (£, e), and solving the cohomological equations as before, we see that we can 
choose go so that s\ becomes independent of x modulo 0((£, e) N+1 ). Repeating this 
argument, we determine successively q 1 ,q 2 , ■ ■ ■ so that 

Sn = s^ + 0^,e) N+1 ), n = 0,1,2,..., 
with = Sn(£,e) independent of x. 
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With Q W = El"o h n q n , we get 

e Q(N> P e e- QiN) = P {N \hD x ,e;h) + R N+1 (x,hD x ,e;h), (3.9) 

where the full symbol of P( N ) is independent of x and the symbol of Rn+i is 
0((h, £, e) N+1 ). Moreover, on the symbol level we have 

P<" +1 >& e; h) - PW(C, 6; fc) = 6, /^ +1 ). (3.10) 

Using this compatibility property we introduce a limiting operator 

P {oo) = P {oo \hD x ,e ] h), (3.11) 

such that 

P<°°\t, e; h) = pW(e, e; h) + 0((£, e, h) N+1 ), 

for all N. Then P (oo) (f , e; h) is well-defined modulo £>((£, e, 

We summarize the discussion of this section in the following proposition. 

Proposition 3.2 Let 

1 



P e ~ p + h Pl + . . . , |f | < 



0(1)' 

Po = Pn(0 + ie{q)(0 + 0(e 2 ) + G(^ +1 ) + cO(^), 
/or some /txed integer N that can be taken arbitrarily large. Then there exist 

N-l 

G^^^^G,-^^), (3.12) 
j'=i 

wift Gj = (9-,(e 2 ) ; 1 < j < iV — 1, ferny homogeneous of degree j + 1 m (f , e) and 
depending analytically on x, and 

N-l 

Q( N \x,Z,e-h) = J2V qj (x,Z,e), (3.13) 
i=o 

where qj are analytic, such that on the operator level, e ad< 9 (iV) e i adG(iV) p e i s of the 
form 

P (N) (hD x , e; h) + R N+1 (x, hD x , e; h), (3.14) 
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where the full symbol of P^ N \hD x ,e; h) is independent of x and 

R N+1 (x^,e;h) = 0((h^,e) N+1 ). 

The leading symbol of P^ N \hD x , e; h) is pf\i, e) = p [N) (£, e) = a ■ £ + le (g)(0) + 
0((£,e) 2 ) = p(£)+ieq(£,e), where 

d^p(0), Re 0^(7(0, 0) are linearly independent. (3.15) 

Remark. Recalling the operator P°°(hD x , e; h) from (|3.1ip . we notice that formally, 
Proposition 3.2 leads to the quasi-eigenvalues 

p(oo) (4 fc -i)-^^) + ^ 00 )' ke7? > ( 3 - 16 ) 

provided that e = 0{h 5 ) and kh = 0{h s ), for some fixed 5 > 0. Here we have also 
used that the space L^(T 2 ), equipped with the L 2 -norm over a fundamental domain 
of T 2 has an orthonormal basis given by 

e k {x) = ei< h ^-^) , keZ 2 . (3.17) 

Remark. The pseudodifferential Birkhoff normal form construction carried out in 
Proposition 3.1 and Proposition 3.2 could also have been done in one stroke, by 
expanding the full symbol of P e in a sum of homogeneous terms in all the variables 
(h, £, e), and making a similar decomposition for the operator G^ N \ We hope to be 
able to return to this idea in a future paper. 

4 The Hamilton-Jacobi equation 

The results of this section will be instrumental in carrying out a complete spectral 
analysis of the operator P £ , in the case when e is sufficiently small but fixed. Here 
we shall do the work that will allow us to construct an additional complex canonical 
transformation which, together with the normal form construction of section 3, 
will allow us to reduce the leading symbol of P e to a function on T*T 2 , which is 
independent of the spatial variables. The construction will be somewhat similar to 
the corresponding constructions in [2E] and |H], and as there, the basic idea is to 
work with cohomological equations of cMype, liberating ourselves from Diophantine 
conditions. 
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Let us recall from Proposition 3.1 that on the principal symbol level, we can 
reduce ourselves to the case of a symbol 

p(x,Z,e)=p( N - 1 \Z,e) + r N (x,Z,e) (4.1) 

in a complex neighborhood of £ = in T*T 2 . Here p^ N_1 '(^, e) is a polynomial of 
degree at most iV — 1 in (£, e), and r^(x, £, e) = e)^)- Assume more generally 

that r = r N = O (e M + \£\ N ^ ioi 2 < M < N . Here M, iV G N can be arbitrarily 

large but fixed. 
Write 

P (N - 1) (t;,e)=p(0+ieq(l;,e), (4.2) 
and recall from Proposition 3.1 that (p is real and that) 

d^p(0), c^Reg(0,0) are linearly independent. (4.3) 

Let ?> be a small parameter and assume that £ = 0(e). We put 

*(£,e)=K0+*«7(£,e)- ( 4 - 4 ) 

We shall try to find a grad-periodic function -^(x) = ij;(x,£) with ^'(a?) = £>(e)j so 
that cp(x,£) = x ■ £ + -^(x) solves the Hamilton- Jacobi equation 

p(x,ip'(x),e)-z({,e) = 0, (4.5) 

so that 

p(£ + + + e) + nvfa £ + e) = z(£, e). (4.6) 

Here we "linearize" the first two terms and get 

(p / (0 + ^(£,e))-«9^ + r^(x,e + ^(x),e) + A(£,^(x),e)^(x)-^(x) = 0, (4.7) 
where 



A{£,V,e)= [ (l-t)(p + ieq)" (£ + tr))dt. 
Jo 



(4.8) 



In order to apply some standard results about non-linear functions of Sobolev 
class functions, we write 

^ = ep, (4.9) 

and (|4.7jl becomes 

(p'(0 + ieq'^; e)) ■ + G(x, f , /4 e, e) = 0, (4.10) 
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where ^ 

G(x,£,ri,e,e) = -r N (x, £ + eq; e) + 7A(£, eq, e)rj ■ rj. (4.11) 
e 

Let yj(x) = yj(x,£; e), j = 1, 2, be linear functions of x 6 R 2 , with 

e ) ■ 9 *Vk = S j,k, ( 4 - 12 ) 
where li = and Z 2 = e). Put 

W * = i ( ej/i + ^ 2 J , (4.13) 



w=- \ey 1 - -2/2 1 , (4.14) 



2 V i 

and 

?/; = , ... , 

2 V i 

so that 

+ e))-&w*=e, (4.15) 

We now look for grad-periodic solutions to (|4.10|) of the form 

p = p peT (x,€,e,e) + b(£,e,e)w*(x,£,£), (4.16) 
where p per is single- valued, and try to find p as a limit p = linij^oo pj, where 

Po = 0, (4.17) 

(^(0 + Hit, e)) ■ d xPj + G (x, £, ( Pj -i)' x , e, e) = 0, j > 1. (4.18) 

We shall estimate the gradient of the solutions to ()4.18|) in the standard if s -norms 
on T 2 for s > 1, and observe that 

ll^ll^~ll(P P er);il^ + |fc|, (4-19) 

uniformly with respect to e, for p of the form ()4.16j) . We shall also use the fact that 
if 

(p't + iecQ ■ d x u = v, 

with u, v periodic, then 

M / I 0(1)11 I 

U> T \H a < \\V Iff.. 
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We notice now that for a given function pj_i of the form ([4.16)1 . ()4.18|) is solvable 
with a solution of the same form, unique up to a constant, and that this equation 
can be decomposed into the following two equations 

eb j + T(G(-^, (Pi-i&ce)) (0) = 0, (4.20) 

and 

(p'(0 + ieq'tiZ, e)) • d x (p peiJ ) + G(x, £, (p^, e,e)-T (G(; £, { Pj „ x )' x , e, e)) (0) = 0. 

(4-21) 

Here we let !Fu(k) stand for the Fourier coefficient at k E Z 2 of the (27rZ) 2 -periodic 
function u. 

In order to treat (|4.18J) with j = 1, we notice that 

|| G(-,£,0,e, e) | = 0(1)= (e M +I N ) . (4.22) 

e 



For pi, we then get 



e M + ~N 



IIPill^ < t ■ (4.23) 

We may remark that for some component of the gradient, we get a better estimate: 

\\p'(0-d x pi\\H°<O(l)- 

For j > 2, we consider in general for u' x , v' x = 0(1) in H s , for some s > 1, 
G(x, £, u' x , e, e) - G(x, f , <, e, e) 
= -{r N (x,£ + eu' x ,e) - r N (x,£ + ev' x ,e)) + [eA(£,er),e)r] ■ r}]"*, =1 + 11, 

6 x 

with the natural definitions of I and II. Writing 

1= / (d ri r N )(x,£ + e(tu x + (l-t)v x ),e)dt •«-<), 



o 



we get 



I || ff .< 0(l)ke M + e N ) ||< -v' x \\ H s 
e 



since by Proposition 2.2 of Chapter 2 in pQ, which extends to the case of vector- 
valued functions, 

\\d n r N {x,Z + e{tu' x + {l -t)v' x ),e) \\ H . < 0(l)l(e M + e N ). 
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Similarly, still assuming that u' x , v' x = 0(1) in H s , we get 

nun*- < o(i)m\u' x \\ H s + \\v' x \\ H s)\\ U ' x -v' x \\ H s. 

Consequently, 

\\G(-,Z,u' x ,e,Z)-G(;Z,v' x ,e,Z)\\ H s (4.24) 
< 0(l)(ke M + e N )+e(\\<\\Hs + \\v' x \\ H s)) \\< - v' x \\ H s. 



e 

We now make the a priori assumption 

|| pf j 11*. = 0(1), j>0, (4.25) 

uniformly in j, e, and notice that it is valid for j = 0,1, by ()4.23|) . if e M ~^ /e, 
~n-i j e < (9(1). Take the difference between the equations ()4.18|) for j + 1 and j: 

(p'(£) + ^(£; e)) • - p,) + (G(x, £, p' p e, e) - G(x, £, p^, e, e)) = 0. (4.26) 

In view of (|4.24|) we get 

II^+i-Pj-IIh* < 0(i) ( — ^~ — + -(II Pj- 11^ + | l^-i ||h«) ) llPj-Pi-i 

(4.27) 

Choose ? with 



< 1, (4.28) 



emin(e, e) 

and notice that this requires that M > 2. Then (J4.27j) implies that 



P'j+i-p'j \\h° < \\\p'j- Pj-i Whs, (4.29) 



2' 

as long as 

|| P ; ||, llp^H^i) 



e M + e N 
ee 



uniformly in j. This holds for j = 1, and as long as we have ([4.29)1 . we can extend 
it to the next j-value. Then (|4.27|) implies that we have a convergent sequence: 



e M + ~N y-l e M + ~N 



p' j+1 ~ p'jWh* < 0(1) . = (4.30) 

1 emm(e, e) J ee 



Summing up, we have 
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Proposition 4.1 Assume that e > 0, e > small, are such that 

< 1. (4.31) 



e M + e N 
emin(e, e) 



Then for £ = 0(e), £/ie Hamilton- J acobi equation ()4.5|) /ias a solution (p(x,£) = 
x ' £ + £p( x ) 0; P °/ ^ e / orm (JHEJ), awd 

||pLH^ = 0((6 M + e^)/(e?)), S >1, 

Essentially, the same iteration scheme as above shows that the solution to (|4.5jl 
is unique up to a constant, if we require that -|| p' x \ \hs -C 1 in addition to ()4.31|) . 



5 Global Grushin problem for e = 0(h 6 ) 

Throughout this section, it will be assumed that e = 0(h 5 ) for some 5 > 0, fixed, 
but arbitrarily small. Our goal is to prove Theorem 1.1 and show that the quasi- 
eigenvalues introduced in ()3.16|) give all the eigenvalues of P e , modulo 0(h°°), in a 
rectangle (|1.25|) . In doing so, we shall only rely on the results of sections 2 and 3. 

Our starting point in this section is the operator P e with a leading symbol (J2.57j) . 
We shall first discuss how to implement the conjugation of P e by the analytic Fourier 
integral operator e ~ lG{N) introduced in Proposition 3.2, which quantizes the com- 
plex canonical transformation k = exp (H G (n)), defined in (|3.4j) . (|3.5j) . In doing so, 
let us consider the IR-manifold, defined in a complex neighborhood of £ = and 
given by exp (H G{N) )(T*T 2 ) C T 2 x C 2 . Here f 2 = T 2 + zR 2 stands for the standard 
complexification of T 2 . It follows from Proposition 3.2 that along this manifold we 
have 

Imx = C(e 2 ), Im£ = C(e 2 ). (5.1) 
Let us introduce the usual FBI-Bargmann transform 

Tu{x) = Ch~ 3/2 J e iip(x ' y)/h u(y) dy, C > 0, 

acting on L^(T 2 ). Here the integration is performed over the whole of R 2 , with 
u being Floquet periodic — see also the discussion in section 3 in [2HI- The phase 
function ip(x, y) = ^(x — y) 2 is defined for x £ T 2 +iR 2 , and the associated canonical 
transformation kt is given by 

T*T 2 3 (y, n) i- (x, = (y - in, v) e T*T 2 . (5.2) 
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The transform k t maps the real phase space T*T 2 to the IR- manifold 

A* : £ = = -Imx, $ (x) = i (Imx) 2 , 

contained inside T*T 2 . Here the zero section in T*T 2 corresponds to Imx = 0. On 
the transform side, the manifold exp (H G (n)) (T*T 2 ) is represented by 

, 2 <9$ e |T , 1 
4 = --^ — , | Imx | < 



i dx ' 1 1 - 0(1)' 

where $ e is a smooth strictly plurisubharmonic function, such that 

V($ e -$ o ) = 0(e 2 ). 

By choosing the undetermined constant in $ e suitably, we can even arrange that 

$ e - $o = 0(e 2 ). (5.3) 

Let now x — xfi mx ) £ Co°> < % < 1, be a standard cutoff function in a 
neighborhood of 0, and consider the IR-manifold given by 

^ ■ * i dx 

where 

$ e (x) = x(Imx)$ e (x) + (1 - x(Imx)) $ (x). 

Here $ e is strictly plurisubharmonic, and if we introduce the globally defined IR- 
manifold 

A e := K?(A it ) C T 2 x C 2 , (5.4) 

then A e agrees with exp (H G ( N ))(T*T 2 ) in a complex neighborhood of the zero section 
and it is equal to T*T 2 further away from this set. Moreover, since along Aj £ , 

we conclude that along A e , we have 

Ime = C(e 2 ), lmx = 0{e 2 ). (5.5) 
Now the Fourier integral operator 

= 0(1) : L 2 (T 2 ) - iJ e (A e ) 
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is such that the action of P t on the space H(A e ), associated to A e by means of the FBI 
transform T, is microlocally near exp (H G (n))(T 2 x {£ = 0}), unitarily equivalent to 
the operator 

e i G(N) P t e~t G(N) : L 2 e (T 2 ) - L 2 e (T 2 ). (5.6) 

From Proposition 3.1 we know that the leading symbol of (J5.6[) is independent of x, 
modulo 0((£,e) N+1 ). 

Let < e'^C 1 be an additional small parameter such that ?S> max(e, h). It is 
then clear, in view of (|3.15|) . that along T*T 2 , in a region where <T< |£| < 1/(9(1), 
we have _ _ 

iRePel^^y or |ImP e - eReF | > (5.7) 

We now notice that the fact that P e= o is selfadjoint in Lg(T 2 ) implies that the 
symbol of 

P - P* 
ImP e = £ % 
2i 

taken in the operator sense on H(A e ), is 0(e) + 0(eh). It follows then from the 
property ()5.5|) of the IR-deformation given by A e that along A e , outside any e- 
neighborhood of the Lagrangian torus exp (Hg n )(T 2 x {0}), the estimates (|5.7|) still 
hold true. 

Corresponding tojdie manifold A t on the torus side, we get a globally defined 
IR-manifold A C T*M, which is e-close to T*M everywhere, agrees with that set 
near infinity, and in a complex neighborhood of Ai, it is obtained by replacing 
exp (ieH G ) o k^ 1 o (k^)-\T*T 2 ) there by 

exp (ieH G ) ° ^ ° (k (7V) ) _1 o exp (H G(N) )(T*T 2 ) = exp (ieH G ) ° 1 o ^A,). (5.8) 

Here we recall that the real analytic canonical transformations Ki and k,( n > have been 
defined in (12. 2|) and (|2.4|l . respectively. (We have also written here exp (ieH G ) for 
the complex canonical transformation identifying A eG and T*M in a neighborhood 
of A x .) " 

Using Propositions 2.3 and 2.5, and taking into account also the conjugation by 
the analytic pseudo-differential operator exp (Q N (x, hD x , e; h)), with Q N defined in 
(|3.13|) . we arrive at the following result. 

Proposition 5.1 Keep all the general assumptions from the introduction, and in 
particular, ()1.13j) . We write F to denote the mean value of q over the Diophantine 
tori Aj, 1 < j < L, and let us make the global dynamical assumption ([1.240 . In 
the case when the H p -flow is completely integrable, we make the assumption (j2.53[) . 
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Recall also the real analytic canonical transformations Kj, 1 < j < L, and k^ n \ 
where Kj is defined in (|2.2j) for j = 1, and is defined in (|2.4|) . Then the 

composed transform 



(iV) o Kj : neigh(A i; T*M) -> neigh(£ = 0, T*T 2 



maps Aj to £ = 0, and /jas t/ie property that when expressed in terms of the coordi- 
nates x and £ on the torus side, the leading symbol p e of P e becomes 



PN,M)+^A^0 + O{e 2 ) + O{i 

with PN,j(0 = a ' £ + 0(i 2 )- Here N > 1 is a fixed positive integer that can be taken 
arbitrarily large. Using these (x,£)- coordinates, we then define for t; small, 



and assume that dpNj(0) = dj = a and Ked(qj)(0) are linearly independent, 1 < 
j < L. Let < 7 C 1 6e such that ? ^> max(e, h). Then there exists a globally 
defined IR-manifold A C T*M and L smooth Lagrangian tori Ai, . . . , Ax, C A, such 
that when p G A is away from any ^-neighborhood of \jfc =x hj in A ; we have 

|ReP e (p)| > — ^-j-t or \lmP e - eReF| > 



0(1) 1 1 ~ 0(1) 

The manifold A is e-close to T*M and agrees with it outside a compact set. We 
have 

P t = 0(1) : H(A, m) -> H{A), 
For each j with 1 < j < L, there exists an elliptic Fourier integral operator 

Uj = 0(1) : H(A) -> ^(T 2 ), 

such that microlocally near Aj, we have 

UP e = (pf\hD x , e; h) + R N+ltj (x, hD x , e; &)) U. 

Here Pj N \hD X} e;h) + RN + ij(x,hD x ,e;h) is defined microlocally near £ = in 
T*T 2 , the full symbol of P^ is independent of x, and 



RN+i,j(x,£,e] h) = 0((h,£,e 



,JV+l^ 
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The leading symbol of Pj N \hD X} e; h) has the form 

a-i + ieF + 0((e, £) 2 ) = P;(£) + e), 

where, as before, pj is real on the real domain, and d^pj(0) and Re d^qj (0,0) are 
linearly independent. We have 

Im pW (f , e; fc) = eRe q 3 (£, e) + 0(e/i) . 

In what follows we shall choose e > e so that 

h l ' 2 - & <e = 0(h s ). (5.9) 

As we shall see, later the lower bound on (Twill have to be strengthened. 
Our goal is to describe the spectrum of P e in a rectangle of the form 



|Rez| < — , 


Im z — eReF 







(5.10) 



when C > is a sufficiently large but fixed constant. To this end, let us introduce 
the quasi-eigenvalues coming from Proposition 5.1 — see also (|3.16j) . 

^U^^^^^-^-^^j^+O^ 1 )), l<j<L, kez\ 

(5.11) 

with h(k—kj/4:) — Sj/2TT = 0(e). In what follows we shall work under the assumption 
that e is bounded from below by an arbitrary but fixed positive power of h, e > h K , 
K ^> 1 is fixed. Let us also remark that it follows from Proposition 5.1 that the 
distance between two neighboring quasi-eigenvalues is > eh/0(l), provided that iV 
is large enough. 

When z G C is in the rectangle ()5.10j) . let us consider the equation 

(P 6 -z)u = v, ueH(A). (5.12) 

At first, we shall derive an a priori bound for the part of u concentrated away from 
the set U^ =1 Aj. 

In what follows we shall write that a function a 6 C°° (A) (also depending on h) 
is in the symbol class SV(1) if uniformly on A, 

V m a = OJr m ), m > 0. 
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We introduce then a smooth partition of unity on the manifold A, 



L 

1 = Yl *3 + + fa,- + i>2,+ + th,- (5-13) 
j=l 

Here < Xj e C{f(A) H * s a cut-off function to an ^-neighborhood of Aj, 

1 < j < L, such that on the operator level we have, 

[P e , Xj ] = 0(h^ 5 ) : H(A) -> #(A). (5.14) 

To construct Xj satisfying (J5.14)) . we use Proposition 5.1 to pass to the torus model 
and take there a symbol of the form Xo (£/?)) where Xo £ C^°(T*T 2 ) is supported 
in a small neighborhood of £ = in J"*T 2 , and Xo = 1 near £ = 0. Conjugating 
Xo{hD x /e) by means of the microlocal inverse of the operator Uj of Proposition 5.1, 
we obtain the cut-off function Xj with the required properties. 

The functions < ipi t ± G «SV(1) in (j5.1rij) are chosen so that ±ReP e > e/0(l) in 
the support of ^,±5 respectively. Finally, we have ± (ImP e — eReF) > ee/0(l) in 
the support of ip 2 ,± £ Co°(A) PI and as in the case of Xji we arrange so that 

in the operator norm, 

A[P e ^ 2 , ± ] = 0(h^ s ), (5.15) 

where A is a microlocal cut-off to a region where |ReP e | < 7/0(1). 
When u satisfies (|5.12|) . we shall prove that 

+ 0(h^ N+1 ^- K - 1 )\\u\\, (5.16) 

where we shall choose N such that (N + 1)5 — K — 1 ^> 1, and provided that the 
lower bound on <T in ()5.9|) is suitably strengthened. Here the norms are taken in 
H(A). When establishing ()5.16|) . we shall first prove that 

< ^^IMI + £ , (^°°)IMI, (5.17) 
e 

which is essentially an elliptic estimate. To prove (|5.17|) . we may follow the method 
of [23]. (See also [T7j.) When M G N, we let 




=: -00 -< ipi -< ■ ■ ■ -< tp 



M 



be a nested sequence of functions in SV(1), supported in a region where ReP e > 
7/0(1). Near the support of ipj we have Re (P e — z) >e/0(l). Using h/e 2 as a new 
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semiclassical parameter that appears naturally on the operator level when reducing 
SV(1) to 5*1(1) by dilation, and applying the sharp Garding inequality we get 

Re((P 8 -z)^u|iM > (^ i j-O(l)^||^i«|| a -O(/» 0O )||tt|r (5-18) 

> ^yll^ll 2 -^ 00 )!^!! 2 , J = 0,...,M. 

provided that h/e 3 <C 1, and we shall even strengthen this further by assuming that 

4 < h s . (5.19) 

Now using that tpj+i) = 0(h°°) in the operator norm, we see that the absolute 

value of the left hand side in ()5.18|) does not exceed 

i; HII ^-u H + ||^ + iw|| 2 + 0(^ oo )||w|| 2 . 

Here we have also used that [P e ,ij)j\ = 0{h/!?). We get 

"^|| 2 < 0(1)|| v (HI tp jU \\ + o(^] \\ip j+1 u\\ 2 + 0(h°°)\\u\\ 2 , 



0(l) llr '~" " - wm-....^-.. ■ - y-p^ 
and therefore, 

ll^ll 2 <^IMI 2 + 0(i^) \\i> j+1 u\\* + o(h°°)\\u\\ 2 . 

Taking into account (j5.19jl . we obtain 

||^H 2 < ^lkl| 2 + 0(^)ll^ + iw|| 2 + 0(^° o )lkl| 2 , 
and combining these estimates for j — 0,1, ... , M, we get 

||^i, + w|| 2 < ^\\v\\ 2 + M (h m )\\iPMu\\ 2 + 0(h°°)\\u\\ 2 . 

The estimate ()5.17j) follows, and by repeating this argument, we get the same esti- 
mate also for ipx-U. We next notice that lmP e = 0(e) on A, and near supp?/>2,+ we 
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have lmP e — eKeF > ee/0(l). An application of the sharp Garding inequality as 
in (|5.18|) gives then for z in (|5.1U|) . 

Im ((P e - z)^ + u\ ^, + u) > 1 1 1 1 2 - 0(h°°) \\u\\ 2 . (5.20) 

The left hand side of (|5.2(J|) does not exceed 

0(1)\\ V || || Ip2,+U\\ + || [P e , Ip2,+]U || || Ip2,+U\\, 

and combining this with ()5.20|) we get 

1 1 1 1 < ^| | v 1 1 + ^| I [P e , fo,+)u 1 1 + 0(h°°) \\u\\. (5.21) 

Now using (|5.15J) together with ()5.17|) and (|5.19J) we get 

1 1 [P„ ^2,+}u 1 1 < 0(h S ) | \V 1 1 + 0(h^ 5 ) 1 1 U 1 1 , 

and we conclude from 1)5.21)1 that 

1 1 ^ u 1 1 < £W | | v 1 1 + ( h (N+l)8-K-l : 



U 



The estimate (J5.10)) follows. 

In the spirit of [23!) as a warm-up exercise, we shall now prove that z in the 
rectangle ([5.10)1 avoids the union of e/i/0(l)-neighborhoods of the quasi-eigenvalues 
()5.1 1|) . then P e — z : H(A,m) — > H(A) is invertible. When doing so, we write, for 
1 < j < L, 

(Pe ~ z)XjU = XjV + [Pe, Xj]u, 

where the if (A)-norm of the commutator term in the right hand side is 0{hS N+ ^ s ) \\u\ 
Applying the operator Uj of Proposition 5.1 we get 

(Pf° + R N+1J - z)U jXj u = U jXj v + T NJ u, (5.22) 

where 

T NJ = 0{h {N+1)S ) : H(A) -> L]{T 2 ). 

Using the fact that Rn+ij(x, £, e; h) = 0((h, e,£,) N+1 ) together with the localization 
properties of Uj X j, we can rewrite ()5.22)) as 

{Pj N) - z)UjXjU = UjXjv + T hNJ u, 
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where T^jvj has the same bound as T^j. Taking an expansion in Fourier series, 

v(k- 0)e k (x), v e Lg, 



V[X) = 

fcez 2 



where e\ are defined in (j3.17|) . we next see that the operator Pj — z acting on 
Lg(T 2 ) is invertible, microlocally near £ = 0, with a microlocal inverse of the norm 
0(1/ eh), provided that z in (|5.1()j) avoids the union of the e/i/(9(l)-neighborhoods 
of the quasi-eigenvalues ([5. lip . It follows that for 1 < j < N, 



XjU\ 



< ^21 1 1 v 1 1 + o^^+ijj-jf-i) \\ u \l e > h K , {N + l)5-K-l^l. 
eh 



Combining this estimate with (|5.16p . we see that the operator P e — z : H(A,m) — > 
H(A) is injective, hence bijective, since general arguments based on the ellipticity 
at infinity (|1.6p . ()1.10|) show that it is a Fredholm operator of index zero. 

We shall now discuss the setup of the global Grushin problem for P e — z, in the 
Hilbert space H(A), which will be well-posed for z varying in the rectangle (|5.1()j) . 
When doing so, we will continue to assume that e > h , when K ^> 1 is fixed. 

With z(j, k), being defined in (|5.1ip . let us introduce for 1 < j < L, 

Mj = # k); \Rez(k)\ < ^, (Im;? - eReF | <gj, 

where C is sufficiently large. Then Mj = O^h" 2 ), and we let k(j, 1), . . . k(j, Mj) G 
Z 2 be the corresponding lattice points, so that 

h (k{j,l) " j) " ^ = 0(c), 1 < / < Mj, 1 < j < L. 

We introduce the operator 

R+ : H(A) -> C Ml x ••• x C Mi , 

given by 

R+u{j){l) = (Uj Xj u\^m), 1<J<L, I <l < Mj. 

Here e&fj-n is as in (|3.17p and the scalar product in the definition of R + is taken in 
L 2 e (T 2 ). Define next 

i?_ : C A/l x ■ ■ ■ x C Ml -> #(A) 
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by 

L Mi 

r-u- = J2Y, u -^^ u i le Hj,D- 

j=l 1=1 

Here Uj 1 is a microlocal inverse of Uj. Let us recall that ey.(x) is micro locally 
concentrated in the region of T*T 2 where £ ~ h{k, — ^) — j-. It follows therefore 
that for I < j < L, 

XjR- =Y, u -WW U i le HM + °( h °°) ■ cMl x • • • x qMl - #( A )- ( 5 - 23 ) 
2=1 

We shall next check that when z G C varies in the rectangle (J5.1U|) . with an 
increased value of C, the Grushin problem 

(P e - z)u + R-U- = v, , g 2 ^ 

has a unique solution («,«_) G H(A,m) x (C A/l x • • • x C A/i ) for every (f,t> + ) G 
H(A) x (C A/l x • • • x C Ml ). Moreover, we have an a priori estimate 

IMI + IKH <^-{\\v\\ + \\v + \\), (5.25) 
ee 

where the norms of u and v are taken in H(A), and the norms of U- and t> + are 
taken in C Ml x ■ ■ • x C A/i . Indeed, we first notice that in view of (j5.16|) and ()5.23|) . 



^ - E « 1 1 < M i + o^^-*- 1 



)(||«|| + ||«_||). (5.26) 



On the other hand, applying Xj an d then Uj, 1 < j < L, to the first equation of 
(JE23J), we obtain 



(Pj* - z)UjXjU + u_(j)(l)e km = UjXjV + 
(^•XjM|e fc(ii0 ) = f+(j)(0, 1 < Z < Mj. 



(5.27) 



Here the Zj(T 2 )-norm of w 3 - is 0(h^ N+1 ^) (|| u|| + ||it_ ||). It follows that for 1 < 
3<L, 

\\Xju\\ + || <^(IMI + IK(j) II) + ^C^(II«II + I|u_||). 
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Here, as in ([5.26)1 . we shall choose N so large that ee > eh > h K+l ^> h^ N+1 ^ 5 , 
which together with ([5.26)1 implies the injectivity, and hence the well-posedness of 
the Grushin problem ([5.24)1 . The bound (|5.25|) follows. 
The solution to (|5.24j) is given by 

u = Ev + E + v + , u_ = E_v + E^ + v + , 

and we recall (see [42j and further references given there) that the eigenvalues of 

P e in the rectangle (|5.1()jl are precisely the values z for which the matrix E |_ G 

£(C Ml x • • • x C Mi , C Ml x • • • x C Ml ) is non-invertible. Arguing precisely as in j2H] 

and j2S] j we find that modulo an error term 0(h^ N+1 ^ 5 ), E |_ is a block- diagonal 

matrix with the blocks E-+(z)(j) G C(C Mi ,C Mi ), 1 < j < L, given by 

E- + (z)(j)(m,n) = {z- z(j,k(j,m)))6 mn , 1 < m < n < Mj. 

The discussion above is summarized in the following theorem. 



Theorem 5.2 Let F stand for the mean value of q along the Diophantine tori Aj, 
1 < j < L. Whenaij andcxij are the fundamental cycles inAj, 1 < j < L, we write 
Sj = (Sij, S2,j) and kj = (k(atij), k(c£2j)) for the actions and the Maslov indices of 
the cycles, respectively. Assume furthermore that e = 0(h s ), 5 > satisfies e > h K , 
for some K fixed but arbitrarily large. Let e > be an additional small parameter 
such that e and 

h l/3-S <? = QQfy 

We next make the global dynamical assumption () 1.24)1 . and assume that the differen- 
tials of the functions Pn,j(0 and Re (qj)(£), defined in Proposition 5.1, are linearly 
independent, when £ = 0, 1 < j < L. In the case when the H p -flow is completely 
integrable, instead of ([1.24)1 we assume ([2.53)1 . Let C > be sufficiently large. Then 
the eigenvalues of P e in the rectangle 



\Rez\ < ^, \lmz - eReF \ < ^ (5.28) 



are given by 



P^^k-^-^e;h^+0(h^% keZ\ 1<J<L. 

Here N G N is such that (N + 1)6 - K - 1 > 1, and 

if%, e; h) = p N>j (0 + ie( qj )(0 + 0(e 2 ) + 0(h). 
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Let us finally recall from ()3.11|) the limiting operator P- ' (hD x , e; h), 1 < j < L, 
well-defined modulo 0((h, e, £)°°) and such that for each iV e N, 



Then it follows from Theorem 5.2 that the eigenvalues of P e in the domain ()5.28j) 
are of the form 



This completes the proof of Theorem 1.1. 

Remark. In the end of this section, we would like to mention that the method of 
Grushin reduction, exploited in this section, has a very long tradition and is closely 
related and essentially equivalent to the so-called Feshbach projection method- 
see [IB] for a recent use of the latter in the context of spectral analysis of the Pauli- 
Fierz Hamiltonian in quantum field theory. We also refer to [12] for a systematic 
presentation of the Grushin method and of some of its applications. In particular, 
in [12], it is explained how the Feshbach method fits into the framework of Grushin 
problems. 

6 Asymptotic expansion of eigenvalues for large e 

In this section we let e be sufficiently small but independent of h and the purpose 
is to compute all the eigenvalues of the operator P e in a domain independent of h, 
thereby proving Theorem 1.2. In doing so, we shall use the results of section 4, and 
we shall also see that a suitable quantum Birkhoff normal form construction will 
allow us to extend the result to cover a certain /i-dependent range of values of e. 

Let us recall from Proposition 3.1 that the leading symbol ()2.57j) of the operator 
()2.56|) can be reduced to the normal form 



with r N (x,£,e) = 0((e,£,) N ), and p N -i(£,e) = p(Q + ieq((,e), with p(f) real, and 
d^p(0), d^R,eq(0,0) linearly independent. Then from Proposition 4.1 we recall that 
the Hamilton- Jacobi equation 



Pj 00 ^, 6; h) = P^'tt, e; h) + 0((h, e, $™). 




P{x, 6, e) = Pn-i{£, e) + r ivO, 6, e), 



(6.1) 



P(x,<p' x ,e) -Piv-i(£,e) = 



(6.2) 



has the solution 



(6.3) 
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for £ = (9(e), e < e, e^/e 2 < 1. Here 1 1 p' x \ \ H , = O^^/e), s > 1. The construction 
in section 4 clearly works the same way for x complex with |Imx| < 1/0(1) and for 
£ complex with £ = 0(e), if we work in if s -spaces on each torus Imx = Const. If 
we normalize the choice of the grad-per iodic function p by putting p(0,£) = 0, then 
p(x,£) becomes holomorphic in x and £, for |Imx| < 1/0(1), £ = 0(e). Write 

<p(z,£,e) = x-£ + i()(x,£,e), (6.4) 

so that _ 

|#r,£,e)| = 0^), (6.5) 

for |Imx| < 1/0(1), £ = 0(e). From the Cauchy inequalities we get 

/~N- 
)";)', ' ' 



3^ = O y-^j (6-6) 

in the same region. 

Consider the actions (independent of x): 

rjj = 1 ^((p(x + 2ire j ,t,e) - <p(x,£,e)) (6.7) 

1 fe N 
= ij + + 2vre„ £, e) - t/j(x, £, e)) = £, + f — 

where e± = (1,0), = (0,1). By the implicit function theorem and the Cauchy 
inequalities, we see that the map £ h- > r)(£), defined for £ = 0(e), is bijective with 
an inverse £(77), and 



r/(£) = £ + 0(^ T J, £(r/) = 77 + j . (6.8) 

Again derivatives of these maps can be estimated by means of the Cauchy inequal- 
ities. Now write ip(x,r),e) instead of ip(x, £,e). We have 

<p(x,ri,e)=x-ri + o(—\ \lmx\ < 1/0(1), \Rex\ = 0(1), |t/| < 0(e), (6.9) 



and ip(x, 77, e) — x ■ 77 is single- valued. 

Consider the canonical transformation 



k = n e : (^</?(x, 77, e), 77) i-> (x, d x ip(x, 77, e)). (6.10) 
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Differentiating the identity ()6.7|) . we get 

d Vk (p(2irej + x, rj, e) = d Vk (p(x, rj, e) + 2n5 jt k , 

so it is clear that n e is a map from a complex neighborhood of the form \lmy\ < 
1/(9(1), rj = 0(e) of the zero section in T*T 2 onto another neighborhood which 
contains a neighborhood of the same form. 

By construction, p o K e (y,T]) = pN-i(£,e), and since £ = 77 + C(e /e), we get, 
expanding the right hand side as a function of (6,77), 

P K e(y, rj) = Pn-i(v, e ) + O (^J^j = P^) + ie( l(v^ t) + {^j^j > ( 6 - n ) 

and this is a function of (77, e) only. 

We next want to implement k by a Fourier integral operator. For that, it 
will be convenient to work in a fixed neighborhood of the zero section. Let P = 
P(x, hD x , e; h) be an /i-pseudodifferential operator with the leading symbol p(x, £, e) 
as in ()6.1|) and assume, as we may, for simplicity, that p(0) = 0, q(0, e) = 0. Put 
h = hfe, so that 

hp(x,hD x ,e;h) = hp(x,ehD x ,e;h). (6.12) 

As an /i-pseudodifferential operator, 7~ l P has a well-defined symbol in a fixed neigh- 
borhood of £ = 0, and the leading symbol will be 

lp(x,e^e) = lp(e£) + lelq^e) + 0(e N ~ l ). (6.13) 

Notice that ? _1 p(e£), e'~ 1 g(?£,e) are uniformly bounded in a fixed complex domain 
when e— > 0, and that 

d^-p(e^), d£-q(e£, 0) are linearly independent for £ = 0, (6.14) 

uniformly when e~— > 0. 
Let 

£/ w (a;) = h~ 2 JJ ei Mx '^ e) - y - v) a(x, 77; dry, (6.15) 

be an elliptic /i-Fourier integral operator, associated to k. The change of variables, 
h = hfe, r] = erj, gives 

Uu(x) = h~ 2 // e^ lip ^ x ' ev '^~ y '^a(x, erj; h)u(y) dydrj. (6.16) 
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From (j6.9|) we get 

1 _ _ fe"-^ 1 

^<p(x,er),e) = x ■ r] + O ( — — 1 , |Imx| < — — , |r/| = 0(1). (6.17) 

The phase is therefore uniformly non-degenerate and the corresponding canonical 
transformation is ^ 

k: (<p' v (x,er),e),r)) (x, ~<f' x (x, efj, e)). (6.18) 

The conjugated operator P = U~ X \PU is a uniformly well-behaved /i-pseudo- 
differential operator with leading symbol obtained from ()6.11|) by division by 7 and 
substitution r\ = erj: 

1 1 _ fe N - l \ 

~p(eff) + ie^q(er], e) + O ( — — 1 =: p(r), e, e), (6.19) 

and this is independent of y. 

We now want to further simplify P by conjugation by an elliptic /i-pseudodiffe- 
rential operator, e A , where A is an /i-pseudodifferential operator of order 0. The 
construction will be uniform in <T with a power degeneration in e, that we shall 
control. First we recall that e A Pe~ A = e adA P = J2 h^a-P- Let the full symbol of 
A be YlkLo h h a k . Then 

e A Pe -A = EEE---E> + ''' +3 " + ' + '(H)--(r d °0 (S) 

fc=0 l=0 ii=0 i fc =0 ' \Al / \A1 / 

oo 

= J^ n sn, 

n=0 

with s = po, si = \H ao p +px = iHp a +pi,... ,s n+1 = iHp a n + s n+1 , where s n+1 
only depends on ao, . . . , a n _i and is a sum of coefficients for h n+1 for terms 

with ji + . . . + jfc + / + A; < n + 1, ji, ... < n. By solving e-degenerated Cauchy- 
Riemann equations, we see that the cij can be successively chosen so that Sj are 
independent of y. 

Assume by induction that Voj = Oie~ x ~ 2 ^\ for j < n — 1 (in a complex domain, 
so that we have the same estimates for the derivatives of Va^). Then the general 
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term in s n+ i is 



\ 2(j 1 +...+j k )+k 

Oil^ 1 -^ 1 . . . e- 1 - 2 ^ = 0(1) i . (6.20) 

Here 

2(ji + ■ • ■ +j k ) + k = 2(j'x + . . . + ]k + k) - k < 2{n + 1 - I) - k = 2n + 2 - 21 - k, 

and so the quantity (|6.20j) is 0{l)e~ 2n except possibly when 21 + k < 2, i.e. when 
= / = or when k = 1, I = 0. In the first case we get the coefficient for /i n+1 in po 
which is 0. In the second case, we get the coefficient for h n+1 in /i J1+1 (iad aji )(po) 
with ji < n, which is 0{l)e' l ~ 2 ^. Here 1 + 2j x < 2n. Thus s n+1 = 0(e~ 2n ), in a 
complex domain. We can choose a n grad-periodic, with iHp a n = — s n+ i+(s n+ i(-, rj)) 
and with Va n = 0(e~ 1_2?1 ). This completes the induction step and we conclude that 
we can find a k with Va k = 0(e~ 1 ~ 2k ) in a fixed complex neighborhood of rj = 
such that if 

Ni-l 

A™ = £ h k a k , 

k=0 

then 

p(m) ;= e AW)^-ii p ^ e -A(^) = (6.21) 

' n=0 

where 

pW) = £>( e ~ 2{n - 1)+ ), (6.22) 

^(^i) _ p(jj^, e) in (j6.19|) and pn ^ = p4°°' > is independent of y and iVi, for n < N\. 
Here iVi can be taken arbitrarily large. 

We shall next look at the quasi-modes. In view of the estimates ()6.22j) on p4 7Vl ' ) , 
we first assume that _ 

-j < 1. (6.23) 
This gives us the quasi-eigenvalues of e^P, 

(S (» - 1) - A ? , e ) + (J^_) , (6 . 24) 

for k E Z 2 with to (A; - ^) - ^ = 0(1). From fEUj) we get the "leading" values 



4 7 27T6 / V V 4 7 27T6 
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and from ()6.19|) we infer that the distance between the neighboring "leading" values 
is > eh/0(l). The ^/-gradient of h n pt\rj, ?, e\ is 0{\)h n e- 2 ^ for 1 < n < N u 
and this is < e, since by (JEH, fc n e -2 ( n-1 > = (^) n e 2 < e 2 . It is therefore clear that 
the distance between neighboring values of 



n=0 ^ ^ ' 



2txT 



, e, e 



is > eh/0(l), and (|6.24|) gives "distinct" quasi-eigenvalues if this minimal separation 
is > ^ iVl+1 /e 2(iVl+1) . We therefore must strengthen (jQED to h^i+i f e 2(N 1+ i) ^ ^ 
i.e.: 



. Ni+l 

h 



< e/i, (6.25) 



or 



h \ Nl+1 eh 



ee 2 I e 



< — . (6.26) 



This will be satisfied if we choose e > h 5 for some fixed (S 6 (0,1/3), e > e, and 
finally N\ large enough. 

Now in ()6.12|) we make the substitution ? — > //?, /i — > h/n, /i ~ 1, and obtain 
the isospectral operator 

-P(x, hD x ) = -P(x,ehD x ) = u— P \ x^ue—D^ 
e e fie \ [i 

with the quasi-eigenvalues 

We deduce as in EE] that for each n = 1, . . . , N%, 



^- n $T ] = $? ) (v,e,e). (6.27) 

We use this to define p4°°^(?7, 1, e) by 

j%°\r],l,e)=?-^(l,Z,ey (6.28) 
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Now recall from Proposition 4.1 that we have the conditions e > e, e N /e 2 <C 1: 

e<?<e 2/Ar . (6.29) 

Then pi 00 ■* (?],?, e) is a well-defined analytic function for \t)\ < 0(1) and by choosing 
to be of the order of magnitude e 2 ^ N , we see that ph°\r], 1, e) is well-defined and 
holomorphic for <^ e 2 ^ N . 
Now restrict e by imposing 

e > ft 1 / 3 " 5 , for some 6 > 0. (6.30) 

Then 

ft ft or 

IP ^ - ^ ^ > 

and the remainder in ()6.24|) can be bounded by 

/ h \ Ni+l 

o(i)( Aj < 

Combining ()6.24|) with (J6.28)) we get that the quasi-eigenvalues of P t take the form 
£ &»p£»> (ft 1, e) + 0(ft 3 ^ +1 )), (6.31) 

n=0 ^ ^ / 71 / 

for fee Z 2 with |fr(Jfe - _ JL| < e 2/W 

Theorem 6.1 Assume that ft 1 / 3-5 < e < eo <C 1, /or some 5 > 0, and let e be an 
additional small parameter such that <T > e and 



~N-3 



e 



2 



< 1. (6.32) 



Recall that F stands for the mean value of q over the invariant tori Aj, 1 < j < L, 
and that Sj G R 2 and kj G Z 2 are the actions and Maslov indices of the funda- 
mental cycles in Aj, 1 < j < L, respectively. In the general case we make the 
assumption () 1.24)1 and when the H p -flow is completely integrable, we assume (|2.53|) . 
Assume finally that the differentials of the functions Pnj(0 o,nd Re (qj)(£), defined 
in Proposition 5.1, are linearly independent when £ = 0. Let C > be large enough. 
Then the eigenvalues of P e in 









Imz — eF\ < — 
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are given by 

*aAO=f> B $?(/i(*-^ fceZ 2 , l<j<L, 

(6.33) 

with 1, e) = 0(e~ 2(n ~ 1 ) + )?~ n ; 1 < j < L, holomorphic for £ = 0(e), and 

^ ) (e,e)=^(0+^(e,e) + O^. 

i/ere pj is real on the real domain, and the differentials of Pj(£) and Regj(£,e) are 
linearly independent when £ = e = 0. The integers N and N x in (16.32(1 and ()6.33() 
can be taken arbitrarily large. 

When proving Theorem 6.1, we assume for simplicity that F = and recall the 
leading symbol of P e from (|2.57j) . 

PN-i{0 + <c<g>(0 + 0{e 2 ) + eO^ 1 ) + O^). (6.34) 

(Here we have replaced A" by N — 1.) Recall also the canonical transformation k 
defined in (13.4(1 and (|3.5|) . and that we have a second canonical transformation given 

by flnm, 

K e : (d v (p(x 1 T],e),rj) -> (x, ^(x, e)) 

defined in a complex domain, with the restriction that |n| = 0(e). Moreover, 
cp(x,T],e) — x ■ r) + 0(e N /e), and we conclude that n e = 1 + (9(? Ar ~ 1 /e), so that 
the differential of K t is equal to 0(e N ~ 2 /e), while the higher order differentials sat- 
isfy V k K e = OiJ N ~ l - k /e), k>2. 

As in section 5, we can represent the real phase space L = T*T 2 by the IR- 
manifold £ = 2 ^fr(x), &o(x) = |(Imx) 2 , with £ = corresponding to Imx = 0, 
and x varying in a neighborhood of {x; Imx = 0} in T 2 + iH 2 . Then L t = K e (Lo) is 
represented by 

t 2d ®t it I /n^ 
£=-—, Imx <0(e), 

z ax 

with $ e - $ , V$ e - V$ = 0(e N - l /e), V fe ($ e - $„) = 0(e N ~ k /e), k > 1. 
Consider the IR-manifold L e represented by 

2<9$ £ . ~ /Imx\ ^ , . / /Imx\\ , . , 

e = = x ( — J $e(x) + ( i - x [—) ) 
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where x £ C^(R?) is a standard cut-off to a neighborhood of 0, so that L e coincides 
with L e near £ = and with L further away from this set. Moreover, 



~N-2 



V($ 6 -$ ), V($ e - $ e ) = O J . (6.35) 

On the other hand, from section 5 we know that |ImP e o k\ along L Q is of the order 
of magnitude eein the region where |£| ~ ?and |ReP e | /e is small. In view of (J6.35)) 
we then have the same conclusion for ImP e o k along L e , provided that 'e N ~ 2 /e <C ee, 
i.e. for 



~N-3 
~1? 



<1, 



which is precisely (|6.32)l . 

The final IR-manifold A e that we choose is then given by 

exp (ieH G ) o (k- 1 ) o (/))-' o k(L £ ) 

in a complex neighborhood of A 1; and we do a similar modification near the other 
tori Aj, 2 < j < L — see also (I5.8J1 . Here we are also using that K(L e ) coincides 
with k(Lo) in the region corresponding to £ ^ suppx (=), £ = 0(e)- Then, along 
A e , |ImP e | > ee/0(l) outside any ^-neighborhood of the union of the A^-'s, 1 < j < 
L, in the energy slice |ReP e | < e/0(l), and moreover the holomorphic canonical 
transformation k o k e maps a region |£| < 0(e) in T*T 2 onto the intersection of A e 
with a similar neighborhood of Ai. 

Using these facts and the normal form (J6.31|) . we get Theorem 6.1 by solving a 
globally well-posed Grushin problem as in section 5. We then take Ni — ► oo, and 
put e = e 1 ^, where N can be arbitrarily large. This completes the proof of Theorem 
1.2. 



7 The dynamical condition and the isoenergetic 
KAM theorem 

7.1 The completely integrable case. 

In the beginning of this section, we shall discuss the case when the Hp-flow is 
completely integrable. To be precise, as in section 2, we shall assume that there 
exists an analytic real valued function / on T*M such that H p f = 0. The sets 
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A a = V 1 (0) n / 1 (a), a & J, then form a singular foliation of p 1 (0), and we get a 
corresponding decomposition of the energy surface, 

P-\0) = \JK- (7-1) 

Here J C R is a compact interval and we know that for each a G J, A a is a 
compact ifp-invariant set, such that we have the continuity property ([2.430 . We 
now introduce a global assumption that for each a G J, except for a finite set 
S C J, the set A a is a finite disjoint union of analytic flow-invariant Lagrangian 
tori, depending analytically on a, A a = U^j"^A a j. Here N(a) is a locally constant 
bounded function on J\S. As we shall see, this assumption is satisfied in the case 
when M is a convex analytic surface of revolution, with the exceptional set S in that 
case corresponding to the endpoints of J, and with N(a) — 1. In order to simplify 
the discussion notationally, in what follows we shall assume that each A a , a G J\S, 
is connected. 

Remark. In some situations, it turns out to be more appropriate to assume that J 
in ([7.1(1 is a circle, or, more generally, a connected graph with finitely many vertices 
and edges, with S then being the set of vertices — see also |TJ],[III|. Assuming that 
J is a graph, the following discussion goes through with minor modifications — see 
also the end of this section. 

It is now classical (see j2J) that locally, near any regular torus A a in T*M, we 
may introduce analytic action-angle coordinates (xi, x% £i, £2) with Xj G S 1 and 
£j G neigh(0,R), j = 1,2, so that when expressed in terms of these coordinates, p 
becomes a function of £1 and £2, P = p(£,i,&)- The rotation number of A a , 00(a), is 
then defined as the ratio of the frequencies 

[%P:M> ( 7 - 2 ) 

viewed as an element of the real projective line. We shall assume that u>(a) depends 
analytically on a G J\S and that a 1— > u(a) is not identically constant on any con- 
nected component of J\S. Following [2 , we recall then that the torus A a C p~ l (0) 
satisfies the isoenergetic condition when the map a ^ uj(a) is a local diffeomorphism. 

When q is a bounded holomorphic function in a tubular neighborhood of T*M, 
that is for simplicity assumed to be real on the real domain, we shall consider the 
limit of 

1 f T 

(q) T = J 1 ex P ( tH p) dt i T ^ 00, (7.3) 

taken along an invariant torus A a , a G J\S. All the results below will be valid if 
(q)t is replaced by (q)r,K, with K G C™. 
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When analyzing (|7.3|) . we switch to the action-angle variables, so that A a becomes 
the standard torus T 2 , H p along A a becomes bid Xl + b 2 d X2 , b\, b 2 G R, and 



q = q(xi,x 2 ). 



Then as in (j^HH , we get 




It follows as in (|2.5()jl that the limit of (q)r along A a is equal to 



q(0)+ E R k > 1 



(7.5) 



0^k,b-k=0 



We conclude that when A a is irrational (in the sense that its rotation number oj(a) = 
b 2 jb\ is irrational), the limit is equal to the average of q over the torus, (q)A a , 
computed with respect to the natural absolutely continuous measure on A a , with 
respect to which the H p -ftow is ergodic. In the rational case, we write u(a) = m/n, 
where m G Z, n G N, are relatively prime, and introduce 



Using the smoothness of q, it follows then from ()7.5|) that the limit of (|7.3|) along a 
rational torus is equal to 



Here the integration measure used in the definition of the torus average is absolutely 
continuous and comes from the diffeomorphism between A a and T 2 , given by the 
action-angle coordinates. 

Let us summarize the discussion above in the following proposition. 

Proposition 7.1 Assume that the H p -flow on T*M is completely integrable and, 
more precisely, that (j7.1|) holds true, with A a being an analytic invariant Lagrangian 
torus when a G J\S, for a finite set S. When q is a real-valued bounded analytic 
function on T*M and a G J\S, we define the closed interval Qoo{o) C R to be the 
range of the limit of {q)x, as T oo, restricted to the invariant torus A a C p~ 1 (0). 
We then have 

• When the rotation number uj (a) ^ Q ; then Qoo{o) = {(<?)a q }- 



k(ui(a)) : 



771 + 77 . 
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Qoo(a) C (q) Aa + O [ ) [-1, 1], k(co(a)) := \m\ + \n\ . 



When u)(a) = m/n, with m G Z, n G N fremg relatively prime, then 

k(u{a)) 
In particular, 

^ |<5oo(a)| < oo. 

a; u(a)£Q 

i/ere |/| stands for the length of an interval /cR. 

When a E S, we also introduce the compact interval Qoo(a) C R, defined as in 
()2.46|) . We shall assume from now on that the function (q)\ a depends analytically 
on a G J\S and that it extends to a continuous function on the whole interval J. 
We shall also assume that a i— ► (g)A a i s n °t identically constant on any connected 
component of J\S. 

Let us next consider the behavior of the interval Qoo(«) corresponding to a 
rational torus A a in a neighborhood of a Diophantine torus. In doing so, we assume 
that A ao , a G J\S, is such that u := u>(a ) is (a, c?)-Diophantine for some a > 
and d > 0, i.e. it satisfies the usual Diophantine condition — cf. (jl.l5|) . 



P 

UJQ 



>^- d , P ez, <? g N. (7.6) 



Let uj = (J (a), a G neigh(ao, R), be rational, and write u = m/n where m G Z and 
n G N are relatively prime. Then 



m m 
= uj = uj — uj + uj . 

n n ' 



and applying ()7.6|) with | = ^, we get 



\UJ — UJq\ > 



n 2+d- 



Therefore, 



k{oj 



l —<0(l)\uj-uj Q \ 1 /^ 



and an application of Proposition 7.1 shows that Q 00(a) is a closed interval such 
that 

Qoo(a) C (q) Aa + O (\uj(a) - u \°°) [-1, 1]. (7.7) 
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This result should be compared with (|2.5ip . 

We shall next consider the set of Lagrangian tori in p -1 (0), whose rotation num- 
bers are uniformly Diophantine and that satisfy a uniform isoenergetic condition. 
Put for a > and d > fixed, 

^ad = a G «/; distfa, 5) > a, |u/(a)| > a, a; (a) > — p G Z, g G N > . 

(7.8) 

In what follows, the dependence on d will not be indicated explicitly, and we shall 
write Q a = fl a ,d- Introduce next 

n a ■.= n a n {a e J\s, |4(g) aJ > «}, (7.9) 

and notice that the measure of the complement of Q a in the interval J, Cf2 a , is 
small, together with a, when d > is kept fixed. We then define the set of good 
values Q a contained inside the closed interval (q)\ a (J), in the following way, 

<?« = C((?>A.(Cn«)). (7.10) 

It follows that the complement of Q a in (q)\ a (J) has a small measure, when a is 
small and d > is kept fixed. We notice also that it follows from the construction 
that when F G Q a , the pre-image 

is a finite set, 

(<?) Aa 1 ( F o) = ... ,a L ; dj eQ, a \ , 
and an application of ()7.7|) shows that 



I a — dj | 



dist (Qoo(a), F ) > ^ , a G neigh(a j; R), 
for any j = 1, . . . L. 

When considering the intervals Qoo(o) for a G J away from the a,-, j — 1, . . . L, 
we notice that an application of Lemma 2.4 shows that the set 

a^noigh({ai,...a i }) 
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is closed. Our final assumption now is that 

sup |<5oo(o)| is sufficiently small depending on a and d. (7-H) 

aeJ 

When a is away from the dj, j = 1, . . . L, using (|7.11jl . we get 

1 



dist (Qoo(a),F ) > 



o(iy 

We summarize the discussion above in the following theorem. 



Theorem 7.2 Assume that the H p -flow on T*M is completely integrable and as- 
sume that the rotation number to •(a) of the invariant tori A a depends analytically on 
a and is not identically constant on any connected component of J\S . When q is 
a real-valued bounded analytic function on T*M , define (q)iv a , a G J\S, to be the 
torus average of q with respect to the natural smooth measure on A a , and assume 
that a (q)A a is an analytic function on J\S, which extends continuously to J and 
which is not identically constant on any connected component of J\S . When a > 0, 
d > 0, let us define next the set Q a C {q)k a (J) according to (|7.10|) . (|7.9J) . Then the 
measure of the complement of Q a in (q)x a {J) is small when a is sufficiently small 
and d is kept fixed. Assume that 

sup |Qoo(a)| 

ae.J 

is small enough depending on a, d. When F e Q a , we have that for any neighbor- 
hood W of the finite set (q)\{Fo) there exists a constant C(W) > such that 

dist (Qoo(a), F ) > ^^y, a e CW. 

Here Qoo{a), a G J, is defined as in (|2.46|) . so that it is equal to the range of the 
limit of {q)r, as T oo, along A a , when a G J\S. 



7.2 Surfaces of revolution. 

We shall now illustrate Theorem 7.2 in the case when M is an analytic surface of 
revolution in R 3 . In doing so, we normalize M so that the X3-axis is its axis of 
revolution, and we parametrize it by the cylinder [0, L\ x S 1 , 

[0,L] x S 1 3 (8,9) i-> (f(s) cos 9 J (s) sin 9, h(s)), 



59 



assuming, as we may, that the parameter s is the arc-length along the meridians, so 
that (/'(s)) 2 + (h'(s)) 2 = 1. A simple calculation then shows that in the coordinates 
(s, 9), the metric on M takes the form 



g = ds 2 + f 2 {s)d9 2 . 



(7.12) 



The functions / and h are assumed to be real analytic on [0,L], and from ^] we 
recall that the regularity of M at the poles is guaranteed by requiring that for each 



and that /'(0) = 1, f'(L) = — 1, which we shall assume from now on. 

Throughout the following discussion, we shall assume furthermore that M is a 
simple surface of revolution, in the sense that < f(s) has precisely one critical point 
s G (0,L), and that this critical point is a non-degenerate maximum, f"(s ) < 0. 
To fix the ideas, we shall also assume that /(so) — 1- F° r future reference, we notice 
that s corresponds to the equatorial geodesic 7e C M given by s = sq, 9 G S 1 . 
This is an elliptic orbit. 

We now come to analyze the geodesic flow on M, which we shall view as the 
Hamilton flow on T*M of the dual form to the metric g. In doing so, we write 



Here a and 9* are the dual variables to s and 9, respectively. Writing out the 
Hamilton equations for the H p — flow, we see next that the functions p and 9* are 
in involution, and we recover the well-known fact that the geodesic flow on M is 
completely integrable. 

When E > and |F| < E l/2 , F ^ 0, it is true that the set defined by 



is an analytic Lagrangian torus sitting inside the energy surface p~ 1 (E). Geometri- 
cally, the torus Ae,f consists of geodesies contained between and intersecting tan- 
gentially the parallels s±(E, F) on M defined by the equation 



k G N, 



/(2*)(0) = f i2k \L) =0, 




(7.13) 



A E , F :p = E, 9* =F, 



f(s±(E,F)) 



F 



E 1 / 2 ' 
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For F = 0, the parallels reduce to the two poles and we obtain a torus consisting 
of a family of meridians. The case \F\ — E 1 / 2 is degenerate and corresponds to the 
equator s = sq, traversed with the two different orientations. 

The principal actions of the torus Ae,f are given by 

h(E,F)= ads = 2 (e-——) ds } (7.14) 

f-'(a) 

and 

I 2 {E,F) = 2ttF. (7.15) 

The functions I±, li depend analytically on E, F ^ 0. 

Let us now restrict the further discussion to the energy surface p _1 (l). We shall 
derive an explicit expression for the rotation number u(a) of the torus A a := Ai )0 , 
7^ a G (—1,1). In doing so, we introduce the action coordinates r\\ = Ii/2tt, 
i]2 = h/^n, and notice that it follows from (|7.14|) and f|7. 15|) that along p _1 (l), r]i 
becomes an analytic function of 772, 

1 r +(a) / rj 2 \ 1/2 

Vi = (find '■= ~ / 1 - Tifr ds > V2 = a, s±(a) := s±(l, a). 

k J s _ { a) V / 2 00/ 

Differentiating next the relation p(ip(r]2), ^2) = 1, we get that the rotation number 
a; (a) is given by —<p'(a), so that 

a r s+(a) 1 / a 2 V 1/2 

y s _ ( a) r(s) V / 2 (s)/ 

We recall that the torus A a C p — 1 (1) satisfies the isoenergetic condition precisely 
when the expression uj'(a), given by 

1 r +(a) j_ ( _^_y 1/2 , f r +(a) j_ ( a 2 y 1/2 , 

(7.17) 

is non-vanishing. 

From now on, we shall assume that the surface of revolution M is such that a; (a) 
in ()7.17|) is not identically constant. As explained in j44j . examples of such surfaces 
are provided by ellipsoids of revolution, with uj'(a) < for oblong ellipsoids, and 
a/ (a) > for oblate ones. The separating case of the sphere is of course degenerate 
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with u'(a) = 0. More generally, we may also remark that the isoenergetic condition 
rules out the case of analytic Zoll surfaces of revolution. 

Let now q be an analytic function on M, which we shall view as a function on 
T*M. We now come to discuss the long time averages of q along the .Hp-flow, (q)r, 
T ^> 1. In doing so, we shall first consider the case when q = qo = qo{s) is a function 
of s only. It follows from ()7.13|) that along an if p -orbit in A a , ^ a G (—1,1), we 
have 

„2 \ 1/2 



m=±2[l -pm> • (7 - 18) 



and then using the rotational invariance of qo, we get 

fT 



(goMm) = ijT q (s(t))dt = (q ) Aa + o(^j , m G A a , T 



oo. 



Here, as we immediately compute from (|7.18)1 . 

(qo)A a = , , , (7.19) 
J (1, a) 

where, in general, for an analytic function if) we write 

J(V, a) = 4>(s) /2 ds, f(s ± (a)) = \a\ . (7.20) 

Js-(a) [J [ S ) - a ) 1 

We remark that (|7.19j) . (|7.20p provide an explicit description of (go)A a which, as 
before, is defined as the average of go over the invariant torus A a . 

The function {qo)\ a depends analytically on a G (—1,1), and it extends to a 
continuous function on the whole interval [—1, 1]. In what follows we shall assume 
that the function a h- > (qo)A a is not identically constant. 

Example. Let 

Qo(s) = (f\s)) 2 = (f"(s )) 2 (s - s ) 2 + 0((s - s ) 3 ), 8 - s . 

Since go is minimal along the equator s = so, it is clear that (qo)A a is not identically 
constant and achieves its minimum for a — ±1. Here we shall compute a leading term 
in the asymptotic expansion of this function, as \a\ —>■ 1. In doing so, following [3] 
and [44] . it will be convenient to make a change of variables on the surface M. We 
introduce 

arcsin/(s), s G (0, s ) 
7r — arcsin f(s), s&(s ,L). 



62 



It follows then that s = c(cosr), where c(x), x G [—1,1], is given by 

(f\ [0 ,s o] )-\VT^), xe [0,1], 



C[X) 



(f\ [soM )-\VT^), [-1,0]. 



Therefore, sinr = f(s), f'(c(cosr))ds = cosrdr, and making the change of variables 
in (JEUD, we get 

r-°H^ sinrfr(cosr) V>(c(x))fe(x) 
J(V>,a)=/ V(c(cosr — ^ ——dr= — - dx. 

Arcsin|a| (sm 2 f - a 2 )V2 J_vT^ (1 - X 2 - Cl 2 ) 1 ^ 

(7.21) 

Here = x/f'(c(x)), x ^ 0, /i(0) = l/|f '(s )| 1/2 ^ 0. 

Taylor expanding h, h(x) = ho + /iix + /i2X 2 + . . ., x — > 0, we get after a simple 
computation, 

J(l, a) = Tth + ^h 2 (l - a 2 ) + 0{{l - a 2 ) 2 ), \a\ -> 1. 

Now 

g (c(x)) = (/'(c(x))) 2 = ^ = lx 2 + 0(x 3 ), 
and using (|7.21j) . we get 

^0,a) = ^-(l-a 2 ) + O((l-a 2 ) 2 ). 

It follows that 

<?o>A a = = i(l - a 2 ) + - a 2 ) 2 ) (7.22) 

is not identically constant. 

Remark. The preceding example is closely related to the computations in Appendix 
Cin 



We shall now introduce a weak angular dependence in the perturbation q, and 
when doing so, for an analytic function qi on M, we set 

q r ,(s,9) = q (s)+T ] q 1 (s,9), < rj < 1. (7.23) 

As before, we are then interested in the limit 

1 f T 

lim - / q v o exp (tH p ) dt, (7.24) 
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taken along the invariant torus A . When analyzing ([7.24)1 . we again switch to 
the action-angle variables, so that A a becomes the standard torus T 2 , H p along A a 
becomes b\d Xl + b 2 d X2 , bj G R, b\ ^ 0, and 

q n (xi, x 2 ) = qi(xi) + r)q 2 {xi, x 2 ). 

It follows from (JZ23) that (fT^Hjl is equal to 

b-k=0 

Let us summarize this discussion in the following proposition, which is just a 
specialization of Proposition 7.1 to the case of surfaces of revolution, with a pertur- 
bation of the form (|7.23|) . 

Proposition 7.3 On the surface of revolution M , let us consider an analytic func- 
tion g^(s,^) = qo(s) + rjqi(s,9), < r] <C 1. When a G (—1,1), we define the closed 
interval Qoo,ri{o) C R to be the range of the limit of (q v )T, as T — > oo, restricted to 
the invariant torus A n Cp _1 (l). We then have 

• When the rotation number u (a) ^ Q ; then Qco,r){a) = {(^Aa}- 

• When uj{a) = m/n, with m G Z, n G N being relatively prime, then 

Qoo,r,{a) C {q v ) Aa + nO ( — ) ] [-1, 1], k(u(a)) := \m\ + \n\ . 



Let us now define the set Q a as in (J7.8)) . with J = [—1, 1] and S = {±1}, and 
define following ()7.9|) . 

^o,« := n a H{ae (-1, 1)\{0}; K(<?o>aJ > «}■ (7.25) 

We then define the //-dependent set of good values Q a%n contained inside the closed 
interval (q v )A a ([-1, 1]), as in (|7.10|) . 

Qa, v = C (((g )A a + v(Qi)A a ) (Cn ,a)) ■ (7-26) 

The complement of Q a ^ has a small measure, when a is small, uniformly in r] small 
enough, and when d > is kept fixed. 

The general assumptions of Theorem 7.2 are verified in this case, and we get the 
following result. 
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Proposition 7.4 Assume that M is a simple analytic surface of revolution for 
which the rotation number defined in (I7.16|) is not identically constant. When 
Qo — %( s ) is an analytic function on M, let us define (qo)A a as in (|7.19|) . (|7.2(J|) . 
and assume that a t— > (qo)A a is not identically constant. Set when < rj -C 1, 

q v (s,6) = q (s) + 77<?i (s,6>), 

and define (q n )i\ a to be the average of q n over the invariant torus A a , with respect to 
the natural smooth measure. When a > 0, let us define the set Q a ^ C (q v )x a ([ — 1> 1]) 
according to (|7.26|) . Then the measure of the complement ofQ a ^ in (q v )A a ([— 1, 1]) is 
small when a and r\ are sufficiently small (and d is kept fixed). When F 6 Q a ri and 
i] is small enough, we have that for any neighborhood W of the finite set (q v )] i 1 (F ) 
there exists a constant C(W) > such that 

dist (Qoo, v (a), F ) > ^|7y, a e CW. 

Here Qoo,r)(a) is the range of the limit of (q^r, as T —> oo, along A a . 

It follows from Proposition 7.4 that for F e Q a ,r)i the assumptions of Theorems 
1.1 and 1.2 (see also Theorems 5.2 and 6.1) are satisfied. An application of, say, 
Theorem 1.1 gives then the complete asymptotic expansions of all the eigenvalues 
of —h 2 A + ieq v in rectangles in the spectral complex plane, that are of the form 

h s eh s 
|Rez-l|<— -r, \lmz - eReFo| < 5 > 0, 

when e = 0(h s ) is bounded from below by a fixed positive power of h. We refrain 
from formulating precisely the final result, as it is immediately obtained from the 
statement of Theorem 1.1. 



7.3 Perturbations of completely integrable systems. 

We now come to discuss the perturbed situation and before doing so, we shall pause 
to recall the statement of the isoenergetic KAM theorem. 
We consider 

PA=Po + Api, 0<A<1, (7.27) 

where po and p\ are holomorphic bounded functions in a tubular neighborhood of 
T*M, that are real on the real domain. We assume that the if po -flow is completely 
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integrable, and let A C p 1 (0) be an analytic //^-invariant Diophantine Lagrangian 
torus. Take then the analytic action-angle transformation 

K : neigh(A , T*M) -»■ neigh(£ = 0, T*T 2 ), (7.28) 

such that k(A ) is the zero section in T*T 2 , and 

p o k ^ = p (0 = « • £ + C(£ 2 ) ; a = a 2 ) G R 2 , 

with a satisfying the Diophantine condition where we assume that Nq > 1. 

Composing p\ in (|7.27j) with the inverse of k in (|7.28j) . we get a function 

pa(z,0 = po(0 + M(^0> 

analytic in a fixed complex neighborhood of £ = 0. An application of the implicit 
function theorem shows that the energy surface p^iO) takes the form £ 2 = 
where / is an analytic function near G R, with /(0) = 0, /'(0) = —ai/a 2 . We 
assume that the isoenergetic condition holds: 

f" (0) ^ 0. (7.29) 

The condition ()7.29|) means that for <C 1, the if po -invariant tori A M : £i = /x, 
£2 = /(a*) can be parametrized by the corresponding rotation numbers /'(/x). 

Theorem 7.5 ("ZTie isoenergetic KAM theorem). Make the isoenergetic assumption 
(|7.29J) and /et us define for a > sufficiently small and d > /ixed, 



i^ a = </i6 neigh(0, R) 
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Assume that < A < « 2 . T/ien there exists a C°°-map 

ty\ : T 2 x neigh(0, R) -> neigh(£ = 0, T*T 2 ) (7.30) 

analytic in the angular variable (uniformly with respect to the other factor) , depend- 
ing analytically on X > 0, and with neigh(0,R) and neigh(£ = 0,T*T 2 ) in ()7.30j) 
uniform in X, such that for each fi G K a , the set 

A M = {^ x (x, /x); x G T 2 } C T*T 2 (7.31) 

zs a uniformly analytic Lagrangian torus C p~^ (0), c/ose to A^ C Pq 1 (0). ITie 
restriction o/exp (tH Px ) to is conjugated to the flow of a constant vector field 
with the rotation number /'(//) on T 2 , fry means of an analytic diffeomorphism, 
depending analytically on X and smoothly on \i G neigh(0,R). Moreover, relative 
to a sufficiently small neighborhood of £ = in T*T 2 , the Liouville measure of the 
complement of the union of the invariant tori A^, /j G f^, m^'fO), is 0(a). 



66 



Remark. The classical theorem of Kolmogorov on the existence of an invariant torus 
for the perturbed Hamiltonian, close to any given Diophantine invariant torus for the 
initial integrable system, is stated and proved in great detail in j3]. The isoenergetic 
version of the Kolmogorov theorem is stated explicitly in but without a proof, 
and a direct argument showing how to derive the isoenergetic version of the theorem 
from the usual one is described in [T7j. The paper j33] gives a complete proof of 
the full KAM theorem and establishes, in particular, the C°°-dependence of the 
KAM tori on the frequencies, in the sense of Whitney. The statement of Theorem 
7.5 above can be extracted from the references [5] (see, in particular, Lemma 2.3 
of that paper), and especially, [H]. The latter paper derives Theorem 7.5 from the 
ordinary KAM theorem of 33j. Finally, we would also like to mention an interesting 
article jH], which provides a complete proof of the isoenergetic KAM theorem. It 
seems, however, that only a continuous, rather than a Whitney-smooth, dependence 
of the KAM tori on the rotation numbers, is obtained in that work. 

Remark. A more refined version of the isoenergetic KAM theorem, which also follows 
from the arguments of jH], gives Cantor families of Diophantine invariant tori with 
prescribed rotation numbers, in energy surfaces p~^ 1 (E), when E G neigh(0, R), with 
an analytic dependence on E. 

We shall apply Theorem 7.5 to the situation described in the beginning of this 
section, with p = p. Let us recall that the energy surface p _1 (0) is foliated by the 
Lagrangian tori A a , a G J\S, and that the rotation number of A a , u(a), is assumed 
to be not identically constant on any connected component of J\S. As before, 
let q be a real-valued bounded analytic function on T*M, and define the function 
(q)A a as the average of q over the if p -invariant tori A a , a G J\S. Assume that the 
assumptions of Theorem 7.2 are satisfied and let F G Q a , with the latter set defined 
in (JD2J), (f7~TU|) . Then 

(^aKFq) = • • • > a L, Oj g n a }, 

and 

dist(Qoo(a), F ) > 1/(9(1), when a is away from the aj, j = 1, . . . , L. (7.32) 

It follows from the equivalence of the assumptions (|1.24}1 and (|2.53jh noticed in 
section 2, together with Lemma 2.2, that for any sufficiently small if p -invariant 
(5-neighborhood W s of Uf =1 A aj in T*M we have, for p G p" 1 ((-5(5), 5(5)))\W S , 
provided that 5(5) > is small enough, 

\(q) S -»i#M-Fo\>^L, iVieN\{0}. (7.33) 
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Let now pi be a bounded analytic and real-valued function on T*M, and let us 
consider 

p x =p + \ Pl , A > 0. 

An application of Theorem 7.5 shows that for A/a 2 small enough, for each 1 < j < L, 
there exists a smooth family of tori A^, a G neigh^-, R), close to A ajJ such that if 
a G O a j, where 



Q a ,j = S a G neigh(aj, R), 



f \ p 
uo(a) 

Q 



a 



then C p x (0) is a Lagrangian torus, on which the H Px -How is quasi-periodic, 
with the rotation number u(a). As before, we may remark that the measure of the 
set neigh(aj, R)\f2 aj - is small when a > is small. 

When 1 < j < L, we define a smooth function of a G neigh ( ctj, R), (q) A u), 

obtained by averaging q over the tori A^. Then we have, in the C 1 -sense, as 
A^O, 



a,A = 



and therefore, when a G neigh (a^, R), we get 



da(q) A u) 

a, A 



a 

>2' 1<J<A 



provided that A is small enough. Following (j7.10J) . let us define 

L 

£ a ,F = CM(?} A ( 3 ) (neigh(a i , R)\fi a j) • (7.34) 

a. A 

Then the measure of the complement of £ a ,F is small with a, and if F £ £ a , 
then for each 1 < j < L, there exists bj G neigh ( a^, R), such that 6 3 - G and 
(q) A U) = F, 1 < j < L. Hence the sets A^ A C J9^ 1 (0), 1 < j < L, are uniformly 

Diophantine Lagrangian tori, and the basic assumptions ()1.13j) - (jl.l5|) are satisfied 
for the f/p A -flow. 

To be able to apply Theorems 1.1 and 1.2 to an operator satisfying the general 
assumptions of the introduction and having the leading symbol p^\ = p + Xpi + ieq + 
C>(e 2 ), we shall now discuss the construction of a weight function G for p 6j \, satisfying 
the conclusion of Proposition 2.3. Outside a small neighborhood of U^ =1 A^\ } we 
shall take G = Gt, p , where Gt iP = Gt is defined in ()1.18j) . using the unperturbed 
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symbol p. When working locally near a fixed if PA -invariant Diophantine torus A^. A , 
which we identify with £ = in T*T 2 , we follow (J2.37)) and set 

G=(1- X ,)Gt, p + X,G. (7.35) 

Here G is defined as in fl2.13j) and < \i <C 1 is to be chosen. An application of 
()2.39j) gives that uniformly in /i > we have 

q - H Px G = (q) + (1 - Xv ) {0{Ti N ) + 0(C° + T~°°)) + <D(T£ N ) + O t (X), 

provided that |T| < On(1) . Here we have also used that 

H Px G T , p = q- (q} T ,K, Px + O t {\). 

As in section 2, we take first N fixed but sufficiently large depending on N\. Choose 
then ji ~ 5, with 5 > sufficiently small but fixed as in ()7.33|) . and put T = fi~ Nl . 
Then in the region where < x M < 1, we have 

q-H Px G-F= (q)(C) - (g)(0) + 0(C N ~ m ) + T (X). (7.36) 

On the other hand, away from the union of the tori A^ A , we are in the region where 
G = Gt, p , and there we have, as A — » 0, 

q - H PX G — F = (q) T>KjP - F + O t {\) - o(l). (7.37) 

Using ()7.36|1 and ()7.37j) . together with ()7.33j) . we conclude that with this choice 
of G, the conclusion of Proposition 2.3 is satisfied for p\, provided that A > is 
taken sufficiently small, depending on the parameters Ni and S. Hence the results 
of Theorem 1.1 and Theorem 1.2 apply in this perturbative situation as they stand 
and give complete spectral results for an operator satisfying the general assumptions 
of the introduction, and with the leading symbol 

Pe,\ = p + \pi + ieq + C(e 2 ), 

in rectangles of the form (jl.25j) and (jl.26|) . for each F G £ a ,F - 

The discussion of this section is summarized in the following theorem. 

Theorem 7.6 Let us keep all the general assumptions on the operator P e from the 
introduction, and write the leading symbol of P e as 

p e = p + ieq + C(e 2 ), 
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near p _1 (0) fl T*M. Let us assume for simplicity that q is real-valued. We assume 
that the H p -flow is completely integrable, so that we have a decomposition (|7.1|) . 
When a £ J\S, we let u(a) stand for the rotation number of the invariant torus A a , 
and assume that a i— ► u(a) is analytic and not identically constant on any connected 
component of J\S. When a £ J, let us also introduce the compact sets Q 00(a) C R, 
defined as in ()2.46|) . and recall from ()1.19|) that 



1 



Im (Spec(P e ) n {z; \Rez\ < 5}) C 



inf \JQ (a) - o(l), sup [J Qoo{a) + o(l) 



aeJ 



aeJ 



as e, h, 8 — > 0. 

We then define a function (q)\ a , a £ J\S, obtained by averaging q over the 
invariant tori A a , and assume that {q)\ a depends analytically on a £ J\S , and 
extends continuously to J. Assume next that a 1— ► (q)x a is not identically constant 
on any connected component of J\S. When a > and d > 0, let T a ,d C U ae jQoo(a) 
be the set obtained by removing from U ae jQoo(a) the following set 



U <2oo(a)]U( U Q^ a] 

dist(a,S)<a / V a£j\S, |w'(a)|<a 



U( U 

1 aGJ\S, |d(g) Att |<a 



U U Qoo(pi) 

\aeJ\s, u(a) is not (a,d)- Diophantine 

Then Theorem 1.1 applies to give complete asymptotic expansions for all the eigen- 
values of P e in a rectangle of the form 



h 6 



h 5 



0(1)' 0(1) _ 



+ ie 



h s „ h s 
,F + 



0(1 



0(1) 



Fn £ T n 



(7.38) 



when 5 > and h K < e = 0(h s ), K 3> 1. analogous statement is obtained 
from Theorem 1.2 in the case when h 1 ^ 3 ^ 6 < e <C 1, S > 0. T7ie measure of the 
complement of T a ^ C U ae j<5oo( a ) small, when a > is small and d is fixed, 
provided that the measure of 

U Qoo(a) ] U(U^( C 
v aeJ\s>( a )eQ / \aeS 



is sufficiently small, depending on a, d. 
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When pi is an analytic function in a tubular neighborhood of T*M, real on 
the real domain, with pi(x,£) = 0(m(Ke(x, £)) in the case when M = R 2 ; and 
Pi(x,£) = C((£) m ) in the manifold case, let P e \, \ > 0, be an operator satisfying 
all the general assumptions of the introduction, and with the leading symbol 

Pe,\=P\ + ieq + (D(e 2 ), p x = p + \ Pl . 

When F G T a ^, let A a . C p _1 (0), a, G J\S, 1 < j < L, be the (a,d)-Diophantine 
tori, with \u'(aj)\ > a, and such that the average of q along A aj is equal to Fq, 1 < 
j < L. When A/a 2 <C 1, we introduce a Cantor family of Diophantine Lagrangian 
tori C p^ 1 (0), close to K aj , 1 < j < L, whose existence is guaranteed by the 
isoenergetic KAM Theorem 7.5. We then define a function a \— > (q) A u) , obtained 

a, A 

by averaging q over the tori A^\, and recall from Theorem 7.5 that (q) A u) extends 

a, A 

to a smooth function on a full neighborhood of aj, 1 < j < L. Let £ a ,F be the set 
obtained by removing from UjL 1 (g) A ( J ) (neigh (oj, R)) the set given by 

a, A 

L 

\^}{q) A U) (a G neigh ( aj, R); u(a) is not (a/2, d) — Diophantine^. (7.39) 
j= i 

Then the measure of the complement of £ ay F is small when a is small and d is 
kept fixed, and when F G £ a ,F an d h K < e = 0{h 5 ), K ^> 1, 5 > 0, Theorem 1.1 
applies and gives complete asymptotic expansions for all the eigenvalues of P tj \ in 
a rectangle of the form (|7.38|) . with F replaced by F, provided that A is sufficiently 
small. A similar result holds by applying Theorem 1.2 to P t) \, in the case when 
h i/3-s < e < e0; < e < 1. 

In the last theorem and elsewhere in the paper, we have tried to find windows in 
the spectral band where all the eigenvalues of P e can be described asymptotically. 
We shall finally discuss a reformulation of Theorem 7.6 which permits us to describe 
a larger subset of the spectral band where the conclusion of the theorem will hold 
uniformly. At the same time, we shall let J be a graph rather than an interval 
(although this is not essential). 

The completely integrable case. We make the same assumptions of complete inte- 
grability as in Theorem 7.2. We also assume that p~ 1 (0) decomposes into a disjoint 
union 

p-\0) = |J A, 
AeJ 
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where A are closed connected if p -invariant sets. We assume that J is a finite con- 
nected graph with S denoting the set of vertices. We assume that the union of edges 
J \ S has a natural real-analytic structure and that every A G J \ S is an analytic 
torus depending analytically on A with respect to that structure. 

We identify each edge of J analytically with a real bounded interval and this 
determines a distance on J in the natural way. As in ()2.43|) . we assume the continuity 
property: 

For every A G J and every e > 0, 3 5 > 0, such that if (7.40) 
A G J, dist (A, A ) < 5, then A C {p G p~ l {0); dist (p, A ) < e}. 

For A G J\S, let uj(A) be the corresponding rotation number, depending analytically 
on A. Assume 

a; (A) is not identically constant on any open edge. (7-41) 

Assume for simplicity that q is real-valued. (In the general case, simply replace q 
below by Reg.) We also assume 

(g)(A) is not identically constant on any open edge. (7.42) 

Here (q)(A) denotes the average of q\.. As before, we also assume that (q) extends 
continuously to J. 

For every A G J, define the closed interval Qoo(A) as in Proposition 7.1, so that 
Qoo(A) = {(q) (A)}, when A G J \ S and uj(A) G" Q. Again, by Lemma 2.4 we know 
that {(A, z); A G K, z G Qoo(A)} is closed whenever K is a closed subset of J. In 
particular, [J AeK Q°° (^) * s com P ac t for every closed subset K C J. As we have seen 
in the case of surfaces of revolution in subsection 7.2, it may happen that 

U ^oo(A) 
Aew- 1 (Q)uS 

has a small measure compared to that of (q){J \ S). 

We shall first define the uniformly good values in R. Fix d > 0. Given a, /3, 7 > 0, 
we say that r G R is (a, j3, 7)-good if the following hold: 

- r is not in the union of all Q 00(A), with dist (A, S) < a or with a; (A) not (a, d)- 
Diophantine (in the sense of (J7.6)0 . 

- If r = (q)(A), then \d A (q)(A)\, \d A uj(A)\ > a. 

- Let (q)~\r) = {A 1; A L }, with Aj G J \ S, dist (A,-, S) > a. Then 

dist(r, |J g oo (A))> 7 . 

AG J, 
dist (A,U^Aj)>/3 
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If r is not (a, ft, 7)-good we say that it is (a, ft, 7)-bad. Choosing successively 
a, ft, 7 sufficiently small, we see that the measure of the set of (a, ft, 7)-bad values 
in R \ UAew- 1 (Q)us < 5°°(^') can ^ e ma d e arbitrarily small. 

T/ie c/ose to completely integrable case. We now replace p by pa as in ()7.27j) and The- 
orems 7.5, 7.6. Near each Aj with (q)(Aj) (a, ft, 7)-good, we apply the isoenergetic 
KAM theorem 7.5, with A small enough depending only on a,/3, 7 and uniformly 
with respect to A, and get a family of KAM-tori A^a^. 
We say that F 6 R is (a, /5, 7, A)-good if: 

- It is a A-perturbation of an (a, ft, 7)-good value r with = {Ai, A^}. 

- There exist KAM-tori A^a close to Aj, that are (a/2, d)-Diophantine, such that 
F = (g) Aj A for j = 1, ...,L. 

Then the set of (a, /?, 7, A)-good values in any bounded interval has a measure 
that tends to that of the set of (a, ft, 7)-good values in the same interval, when A 
tends to zero. Moreover the conclusions of Theorems 1.1 and 1.2 hold uniformly 
when A is small enough depending on a,ft,j, and F is (a,ft,j, A)-good. 

8 Barrier top resonances: the non-resonant case 

The purpose of this section is to illustrate Theorem 1.2 by applying it to the problem 
of distribution of semiclassical resonances for the Schrodinger operator. The discus- 
sion here will be analogous to the corresponding treatments in j2H], EH] and [2~^j . 
and for this reason, the following presentation will be somewhat less detailed. 
Consider 

P = -h 2 A + V(x), P (x,0=e + V(x), (x,OeT*R 2 , (8.1) 

where V is an analytic potential satisfying satisfying the general assumptions of 
section 7 of [2S1, that allow us to define the resonances of P in a fixed sector in the 
fourth quadrant. Let V(0) = E , V'(0) = and assume that we are in the barrier 
top case, so that V"(0) < 0. In this section, we shall also assume for simplicity 
that V is an even function. The Taylor expansion of p in a suitable system of linear 
symplectic coordinates then takes the form 

p(x, Z)-E = fl y (3 - *?) + P4(*) + • • • , (8-2) 
i=i 

where Aj > 0, j = 1, 2, and p^ is a homogeneous polynomial of degree 4. 
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We shall assume that we are in the non-resonant case so that 

A ■ k ^ 0, O^keZ 2 . (8.3) 

In this case, the result of j2E] gives all resonances of P in the open disc D(E , h s ) = 
{z G C; \z — Eq\ < h s }, for any 5 > 0. They are given by 

z k = E + f(h(k-e )h;h), k G N 2 , 

where 6q G is fixed, and f(8; h) is a smooth function of 6 G neig h(0,R 2 ), with 

/(0; fc) ~ f (6) + hh{6) + .... We have 

f (9) = -i\-6 + O{6 2 ). 

The purpose of this section is to show how Theorem 1.2 allows us to go further away 
from the real axis and obtain a description of resonances that are at a distance ~ eo, 
< 6q <C 1, away from the real axis. 

As in |24j . in order to study the resonances of P near i? , we perform the 
complex scaling, which near (0, 0) is given by x = e t7T ^x, £ = e _l?r//4 £, x, £ G R 2 . 
Using (|8.2|) and dropping the tildes from the notation, we get a new operator with 
the leading symbol 

- (P2(x, - zp 4 (z) + ...)=: -q(x, 0, (x, - (0, 0). (8.4) 

Here 

P2(x,£) = Ey(^ 2 + £J) ( 8 - 5 ) 

is the harmonic oscillator. We shall be interested in eigenvalues E of the operator 

Q{x,hD x ;h) = q{x,hD x ) + O(h), 

with \E\ ~ e 2 , < f < 1. An additional rescaling x = ey, h 5 < e < 1, < 5 < 1/2, 
gives 

\Q{x, hD x ; h) = \Q(e{y,hD y ); h), h = ^<£l. 
Viewed as an /i-pseudodifferential operator, e~ 2 Q(x, hD x ; h) has the leading symbol 

p 2 (x,0-^ 2 P^) + O(e 4 ), 
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to be considered for (x,£) in a bounded region. 

The if P2 -flow is completely integrable, and as in (2B], we introduce the action- 
angle coordinates Ij > and Tj e R/27rZ, j = 1,2, for p 2 , given by 

Xj = y2Lj cos Tj , £j = — y / 2ij sin Tj . 

We also remark that the non-resonant condition (J8.3|) implies that the H p2 -&ow is 
ergodic on each invariant torus A M C p 2 given by 

Therefore, the flow average 



i r 

(p 4 )t = ^ y p 4 o exp (t/Tp,,) dt 



along the torus A M converges to the torus average of p±, (p^) = (^4)00, as T — ► 00. 

When computing the flow and torus averages, it will be convenient to work in 
the symplectic coordinates (y, 77) given by 

In these coordinates p 2 = Ylj^^jUjVjy an d the flow is given by 

exp (tH P2 )(y,r)) = (e itAl yi , e<% 2 , e^ 1 ^, e -'% 2 ). 

Writing 

^ Q = a ^y k V a ~\ l«l = 4, a fcQ = — — ( Jfe ) ' 
we immediately see that as T — ► 00, the flow average (x a )x converges to 



0<k<a 
2k=a 



It follows that if 

Pa{x) = ^2 Va%C 

|a|=4 
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then 

(P4)(y,V) = — (6^(4,0)1/^1 + 4 ^(2,2)2/l2/2^2 + 6^(0,4)2/2^) • 

This function is in involution with p 2 , and can be expressed in terms of the action 
variables as follows, 

(PA) = T (6^(4,0)A 2 + 4<y {2,2)hh + 6^(0,4)^1) • 

Now P2 = Xih + A2/2, and therefore we see that for /1/2 7^ 0, the differen- 
tials dp 2 = (Ai,A 2 ) and d(p 4 ) = -(3f( 4>0 )/i + V( 2 ,2)h, 3w( ,4)-^2 + V(2,2)h) are linearly- 
dependent precisely when 

(Al^(2,2) - 3A 2 f(4,0))A = (A 2 ^(2,2) - 3AlU( 0) 4))i2- (8.6) 

When I\ = or I 2 = 0, the question of linear independence of the differentials should 
be examined directly in the (y, -^-coordinates, and we then easily see that we have 
the linear dependence when /1/2 = 0. We also compute that (p^) restricted to p^ 1 ^) 
has the critical values A\ = fi^o^A^ 2 , corresponding to I\ = 0, A 2 = ft^ojA^ 2 , 
corresponding to I2 = 0, and a third value A 3 which occurs when the line (J8.6|) 
intersects the quadrant l\ > 0, I 2 > 0, i.e., when 

(Al^(2,2) - 3A 2 ^(2,2))(A2^(2,2) ~ 3AiU( ,4)) > 0. 

We are now in a position to apply Theorem 1.2 to the operator 

\q(x, hD x ) - 1, 

with e 2 considered as a small perturbation parameter and using h = 4r as the new 
semiclassical parameter. Here we should also remark that the Birkhoff normal form 
construction of section 3 goes through in the present case assuming only (|8.3|) and 
that the Diophantine condition is not required — see also [37] for the classical Birkhoff 
construction in the non-resonant case, near a stable equilibrium point. 

Theorem 8.1 Consider the operator P in (jO| with the leading symbol 

P (x, o - e = y (£f - +4 + o(x e ), (x, o - (o, o). 
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Assume that the non-resonance condition (A, k) ^ 0, 7^ k G Z 2 , holds. Let us put 
Ai = 0, A 2 = |Aj" 2 . Assume that < S < \. Then the resonances of P in the 
domain 

2 

{z G C; h s < \z - E \ < e }\ [j{z G C; \~Rez - E - Aj \lmz\ 2 \ < r] |Im^| 2 }, 

i=i 

(8-7) 

where < eo = eo( r ?) 1 and t] > is arbitrary small but fixed, are given by 

n=0 ^ ' 

with p^°\^, e) analytic in £ G neigh(0, C 2 ) and smooth in e G neigh(0, R). VFe aai>e 
5 G R 2 and a G Z 2 are /ized, and we choose e > wif/i ji? — i?o| ~ e 2 . 

Remark. Combining Theorem 8.1 with the arguments of section 7 of [2E], we obtain 
that the result of j2H] extends to a set of the form 

2 

D(E , e )\\J{z-Rez-E Q e -[A, - r), A, + ??](Imz) 2 } . 
3=1 

Here D(E Q , eo) = {2, \z — E \ < eo}, with < eo = eo{v) ^ 1> an d 77 > is arbitrary 
but fixed. 
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